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1.  INTRODUCTION 


1.1  OVERVIEW 

A strapdown  micro-navigator  (MICRON)  is  presently 
under  development  by  AFAL/Autonetics  for  applications  requir- 
ing a small,  reliable,  moderately  accurate  inertial  navigator. 
The  entire  set  of  inertial  sensors  in  MICRON  consists  of  two 
Micro-Electrostatic  Gyro  Accelerometers  (MESGA's).  As  part 
of  the  development  of  MICRON,  a strapdown  inertial  navigator 
designated  as  N57A  is  being  fabricated  and  tested  (Ref.  6). 
This  system  uses  three  accelerometers  for  measuring  specific 
force  and  two  Micro-Electrostatic  Gyros  (MESG's)  for  maintain- 
ing an  inertial  reference.  The  physical  characteristics  of 
the  MESG  are  identical  to  those  of  the  MESGA.  The  N57A  is 
intended  to  demonstrate  MICRON'S  function  and  performance, 
but  not  MICRON'S  ultimate  small  size. 

TASC  has  been  performing  analytic  studies  in  support 
of  the  MICRON  development  program  since  March  1971  through  a 
sequence  of  contracts  with  AFAL.  Previous  investigations 
have  looked  chiefly  at  the  MESG  instrument,  developing  mod- 
els to  describe  the  MESG  drift  rate  and  attitude  readout  (ARO) 
errors  (Refs.  1,  2,  and  3)  and  evaluating  UNICAL,  a procedure 
proposed  by  TASC  to  simultaneously  calibrate  both  the  drift 
rate  and  ARO  error  models  (Refs.  3 and  4).  This  report 
documents  a set  of  system  level  analyses  that  have  been  con- 
ducted to  support  the  development  of  MICRON.  These  analyses 
were  conducted  in  three  areas : direct  simulation  of  MICRON 

system  performance  using  the  MICRON  Computer  Simulation  Pro- 
gram (MCSP) , derivation  of  error  analysis  equations  and 
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development  of  a computer  program  implementing  these  to  con- 
duct  statistical  performance  studies,  and  support  of  the 
N57A  test  program. 

The  MCSP , which  simulates  directly  the  mechanization 
equations  that  are  implemented  in  MICRON  (see  Fig.  1.1-1), 
was  developed  during  previous  MICRON  studies  and  is  described 
in  detail  in  Ref.  3.  The  MCSP  was  used  to  investigate  the 
navigation  performance  of  the  MICRON  system  during  four  dif- 
ferent cases: 

• Keducea  preload  charge  levels  on  the  MESG. 

• Severe  vibration  environment. 

• High  vehicle  roll  rates  and  accelerations 

• Instrument  errors  due  solely  to  miscali- 
bration  (i.e.,  assuming  no  errors  due  to 
turn-on-to-turn-on  repeatability ) . 

Plots  of  the  navigation  errors  occurring  in  each  of  these 
cases  were  generated. 


performance 

Figure  1.1-1  MCSP  Flow  Diagram 
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To  complement  the  MCSP  and  to  enhance  the  total  cap- 
ability to  simulate  MICRON  performance,  a covariance  error 
analysis  computer  program  was  developed.  The  covariance 
error  analysis  technique  allows  the  generation  of  the  sta- 
tistics describing  MICRON'S  navigation  performance  for  an 
ensemble  of  trials  over  the  same  trajectory  in  a single 
computer  run.  (To  conduct  statistical  performance  studies 
with  the  deterministic  MCSP  would  require  costly  Monte  Carlo 
analysis.)  This  report  derives  the  MICRON  error  equations, 
which  model  the  propagation  of  instrument  errors  into  sys- 
tem navigation  errors,  and  describes  the  implementation  of 
these  equations  in  the  error  analysis  computer  program. 

The  results  of  a study  to  investigate  the  error  propagation 
characteristics  of  those  error  sources  that  are  unique  to 
MICRON  are  also  presented. 

The  final  activity  discussed  in  this  document  is 
the  TASC  support  of  the  N57A  test  program.  The  Test  Plan 
developed  by  Autonetics  (Ref.  20)  is  reviewed  and  recom- 
mended improvements  are  suggested.  The  results  of  the 
laboratory  and  van  tests  on  the  N57A  system  are  evaluated 
and  the  system  performance  is  compared  with  the  specified 
performance  values  (Ref.  5). 


1.2  OUTLINE  OF  THE  REPORT 

Chapter  2 presents  the  results  of  the  MCSP  studies. 
Section  2.1  briefly  describes  the  MCSP  while  Sections  2.2, 
2.3,  and  2.4  present  the  simulated  MICRON  navigation  per- 
formance for  the  four  cases  that  were  investigated.  A brief 
overview  of  the  mechanization  equations  is  presented  in  Sec- 
tion 3.1  while  the  error  equations  which  describe  the  MICRON 
system  are  derived  in  Section  3.2.  Section  3.3  discusses  the 
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error  analysis  computer  program  and  presents  the  results  that 
were  generated  with  it.  The  N57A  Test  Plan  is  reviewed  in 
Section  4.1,  and  the  results  of  the  test  program  are  discussed 


in  Section  4.2. 
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The  MICRON  Computer  Simulation  Program  (MCSP),  de- 


scribed in  Ref.  3,  was  used  to  investigate  the  effects  of 
severe  environments  and  of  several  proposed  MESG  variations 
on  MICRON  navigation  performance.  The  environmental  stud- 
les  consisted  of  determining  the  performance  degradation  in- 
duced by  vibration  inputs  and  by  high  vehicle  roll  rates. 

The  MESG  variation  studies  evaluated  the  predicted  naviga- 
tion  performance  for  cases  when  the  gyro  preload  charge 
level  was  reduced  and  when  there  were  no  gyro  repeatability 
errors  (i.e.,  the  only  MESG  errors  were  due  to  the  calibra- 


tion procedure). 

A brief  description  of  the  MCSP  is  presented  in  Sec- 
tion 2.1.  The  results  of  the  MESG  variation  studies  are 
given  in  Section  2.2,  while  the  vibration  and  high  roll 
rate  results  are  presented  in  Sections  2.3  and  2.4,  re- 
spectively . 

2.1  MCSP  DESCRIPTION 

The  MICRON  Computer  Simulation  Program  (MCSP)  is  a 
flexible  research  tool  for  investigating  the  performance  of 
the  MICRON  system.  The  program  is  a direct,  deterministic 
simulation  which  allows  study  of  the  parametric  behavior 
of  the  MICRON  system,  its  hardware  and  its  software,  under 
a variety  of  dynamic  conditions.  The  MCSP  is  a direct  sim- 
ulation since  it  duplicates  the  mechanization  equations 
used  in  the  N57A  (a  developmental  system,  using  three 
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accelerometers  to  measure  specific  force  and  two  MESG's  to 
maintain  an  attitude  reference,  which  demonstrated  MICRON'S 
function  and  performance,  but  not  its  ultimate  small  size). 

It  is  a deterministic  simulation  since  the  program  input 
parameters  describing  the  system  error  sources  are  "sample 
inputs;"  and,  the  simulation  results  generated  by  the  MCSP 
are  the  "sample  outputs"  corresponding  to  these  inputs.  To 
perform  a statistical  analysis  of  predicted  MICRON  naviga- 
tion performance  using  the  MCSP,  it  would  be  necessary  to 
make  a large  number  of  simulation  runs  drawing  the  input 
parameters  from  appropriate  distributions  and  average  the 
entire  ensemble  of  results  that  were  generated  (i.e.,  to 
perform  Monte  Carlo  analysis). 

The  MCSP,  diagrammed  in  Fig.  1.1-1,  has  been  devel- 
oped in  a modular  form  to  facilitate  modifications  and  addi- 
tions as  the  MICRON  system  design  is  altered  or  new  areas 
require  investigation.  The  principal  program  modules  include 
a trajectory  generator/truth  model  (TG/TM),  several  subrou- 
tines which  model  the  inertial  instruments,  and  a module  that 
duplicates  the  system  mechanization  equations.  The  TG/TM 
supplies  a dynamic  motion  environment  with  which  to  exercise 
the  MCSP  and  includes  an  accurate  model  of  position  and  vel- 
ocity against  which  to  evaluate  the  MICRON  computation  of 
these  quantities.  The  TG/TM  allows  the  selection  of  several 
different  trajectories  (see  Ref.  3)  for  use  in  MICRON  evalua- 
tion . 


Incorporated  into  the  MCSP  is  a detailed  simulation 
of  the  Micro-Electrostatic  Gyro  (MESG)  which  includes  the 
prominent  drift  rate  and  attitude  readout  (ARO)  error  mech- 
anisms. As  described  in  detail  in  Ref.  3,  only  residual 
(those  which  remain  following  compensation)  MESG  errors  are 
used  in  the  MCSP.  This  procedure  simplifies  the  simulation 
program  by  eliminating  the  necessity  of  implementing  the  full 
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MESG  compensation  formulas.  Parameters  descriptive  of  the 
MESG  residual  errors  for  the  case  to  be  studied  are  part  of 
the  MCSP  input  data.  For  the  current  studies,  the  accel- 
erometers have  been  modeled  as  ideal  instruments. 

The  execution  rate  of  the  MCSP  is  set  by  three  pa- 
rameters: fast  loop  frequency  (FLF),  slow  loop  frequency 

(SLF) , and  truth  loop  frequency  (TLF) . The  first  two  are 
MICRON  system  parameters  which  indicate  how  often  the  iner- 
tial sensors  are  sampled  and  how  often  the  navigation  equa- 
tions are  solved.  For  the  MICRON  system  these  are  64  and 
16  times  per  second,  respectively.  The  truth  loop  frequency 
indicates  how  often  the  navigation  reference  quantities  (i.e., 
true  position,  velocity,  and  attitude)  are  updated  by  the 
trajectory  generator.  All  three  iteration  rates  may  be 
varied  independently  in  the  MCSP.  This  desirable  feature 
provides  the  capability  to  investigate  different  solution 
rates  for  the  mechanization  equations  and  to  vary  the  truth 
loop  frequency  as  a function  of  the  trajectory  allowing  ef- 
ficient , yet  accurate,  reference  computations  for  dynamic 
environments  of  varying  frequency  content . 

Because  it  duplicates  the  MICRON  mechanization  and 
consequently  requires  a large  amount  of  computer  time  to 
execute,  the  MCSP  is  most  useful  for  studies  which  can  be 
conducted  using  short  duration  trajectories  and  which  re- 
quire the  level  of  detail  which  the  MCSP  offers.  Most  long- 
term navigation  performance  investigations  do  not  require 
the  detail  of  a direct  simulation.  For  these  studies,  a 
covariance  error  analysis  program,  which  allows  a statisti- 
cal characterization  of  the  system  performance,  is  a more 
appropriate  analytic  tool.  The  preliminary  development  of 
a covariance  error  analysis  program  for  MICRON  is  described 
in  Chapter  3. 
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2.2  MESG  VARIATION  STUDIES 


The  MESG  variation  studies  investigated  the  effect 
on  MICRON  navigation  accuracy  of  four  different  levels  of 
MESG  errors.  These  residual  error  levels  represent  the  er- 
rors for  various  operational  or  design  variations  potentially 
available  to  the  MICRON  system  designers.  For  each  case,  the 
MCSP  input  parameters  that  model  the  MESG  were  varied  to  yield 
the  appropriate  residual  drift  rate  and  ARO  errors,  and  the 
simulation  program  was  exercised  to  determine  the  resulting 
navigation  errors.  Table  2.2-1  presents  the  residual  values 
that  were  used  in  these  studies  to  represent : 


Perfect  MESG's  (i.e. , no  drift  rate 
or  ARO  error). 

Nominal  MESG's 

MESG's  with  errors  due  only  to  mis- 
calibration 

MESG's  with  the  preload  charge  level 
reduced  by  30% 


The  values  chosen  to  describe  the  drift  rate  and  ARO  error 
residuals  were  based  on  the  MICRON  performance  specifications 
and  on  data  collected  by  Autonetics  during  MESG  testing. 


TABLE  2.2-1 

RESIDUAL  ERRORS  FOR  MESG  VARIATION  STUDIES 


MESG  VARIATIONS 


Perfect  MESG's 

Nominal  MESG's 

MESG's  with  only 
calibration  errors 

Reduced  Preload  MESG' 


RESIDUAL  DRIFT  RATE  RESIDUAL  ARO  ERROR 
(deg/hr)  (mrad) 


0.0100 


0.0092 


0 . 0070 


0.100 


0.078 


0.110 


The  residual  error  levels  selected  for  the  nominal 
MESG  and  for  the  MESG  which  exhibited  calibration  errors  only 
are  the  targets  called  out  in  the  MICRON  specification  (Ref.  5) 
Since  Autonetics  has  demonstrated  (Ref.  6)  that  these  target 
MESG  error  residuals  are  attainable,  their  use  in  this  per- 
formance evaluation  study  is  reasonable.  MESG  reduced  pre- 
load studies  conducted  by  Autonetics  (Refs.  7 and  8)  indicated 
that  using  a preload  charge  level  equal  to  70%  of  the  nominal 
value  leads  to  a 30%  reduction  in  drift  rate  residuals  and  a 
10%  increase  in  residual  ARO  errors. 


Figure  2.2-1  presents  the  trajectory  which  was  used 
for  the  MESG  variation  studies.  The  aircraft  speed  for  this 
5 min,  constant  altitude  (24,000  ft)  trajectory  is  600  ft/sec 
During  the  turn,  the  heading  rate  is  0.75  deg/sec.  Since  the 
MICRON  alignment  procedure  was  not  investigated  in  these  stud 
ies,  the  navigation  initial  position,  velocity  and  attitude 
errors  were  set  to  zero.  A more  detailed  description  of  this 
trajectory  is  presented  with  the  MCSP  development  in  Ref.  3. 


LEVEL 

FLIGHT 


FLIGHT 


Figure  2.2-1  Trajectory  for  MESG  Variation  Studies 


The  time  histories  of  the  MICRON  navigation  errors 
for  the  MCSP  run  using  nominal  MESG  error  levels  are  pre- 
sented in  Fig.  2.2-2.  Examination  of  these  results  indi- 
cates that  the  turn  in  the  trajectory  induces  a MICRON 
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Single  Trial  Monte  Carlo 
MICRON  System  Performance 
With  Nominal  MESG  Errors 


Figure  2.2-2 
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velocity  error.  This  error  results  from  misresolut ion  of 
sensed  acceleration  due  to  ARO  errors  and  misalignments  be- 
tween the  true  navigation  frame  and  the  computational  frame 
caused  by  drift  of  the  MESG’s  (see  Ref.  9 lor  a discussion 
of  error  propagation  in  inertial  navigation  systems).  At 
the  completion  of  the  turn  when  the  aircraft  acceleration  is 
zero,  the  velocity  error  behavior  will  exhibit  the  classi- 
cal 84-min  period  Schuler  oscillation,  although  the  trajec- 
tory is  too  short  to  clearly  demonstrate  this  behavior. 

The  navigation  errors  for  each  of  the  four  cases  in- 
vestigated in  the  MESG  variation  study  have  similar  trajec- 
tory-induced time  histories  to  those  presented  in  Fig.  2.2-2 
for  the  nominal  MESG  errors.  The  magnitudes  of  these  errors 
vary  indicating  the  navigation  performance  attainable  with 
each  of  the  four  levels  of  MESG  errors.  Table  2.2-2  pre- 
sents the  terminal  navigation  errors  for  the  5-min  trajec- 
tory diagrammed  in  Fig.  2.2-1.  The  errors  for  the  perfect 
MESG  run  indicate  the  navigation  errors  due  solely  to  mech- 
anization of  the  MICRON  navigation  equations  (iteration 
rates,  integration  algorithms,  etc.).  Table  2.2-2  shows 
that  mechanization  errors  do  not  dominate  MICRON  perform- 
ance for  this  trajectory  since  they  account  for  less  than 
10%  of  the  navigation  error  when  nominal  error  levels  are 
used  to  describe  the  MESG. 

It  is  important  to  note  that  the  results  presented 
in  Table  2.2-2  are  single  trial  Monte  Carlo  runs  (one  set  of 
sample  inputs  producing  a single  set  of  sample  outputs)  and 
as  with  single  samples  of  a random  process,  care  must  be 
exercised  when  interpreting  the  results.  (Table  2.2-2  does 
not  present  statistical  performance  predictions  for  the  pro- 
posed MICRON  configurations.)  In  generating  the  results 
presented  in  Table  2.2-2,  it  was  assumed  that  the  trajectory 
was  preceded  by  a perfect  alignment;  and  consequently,  the 
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TABLE  2.2-2 

SINGLE  TRIAL  MONTE  CARLO  TERMINAL  NAVIGATION  ERRORS  FOR 

MESG  VARIATION  STUDY 


MESG 

VARIATION 


NORTH  EAST  RADIAL  NORTH  EAST  RADIAL 

POSITION  POSITION  POSITION  VELOCITY  VELOCITY  VELOCITY 

ERROR  ERROR  ERROR  ERROR  ERROR  ERROR 

(**)  (ft)  (ft)  (ft/sec)  (ft/sec)  (ft/sec) 


Reduced  Preload 
MESG ' a 


Perfect  MESG's 

m 

- 3.8 

mm 

-0.010 

-0.034 

Nominal  MESG *8 

20.8 

-60.3 

63.8 

-0.226 

-0.525 

MESG 1 s with 
calibration 
errors  only 

14.4 

-54.1 

56.0 

-0.228 

-0.412 

-0.105  -0.477  0.488 


terminal  errors  are  smaller  than  would  be  expected  from  the 
actual  MICRON  system  since  the  propagation  of  initial  align- 
ment errors  into  system  performance  errors  has  not  been  in- 
cluded. For  these  reasons,  the  conclusions  developed  from  the 
results  presented  in  Table  2.2-2  will  be  based  on  the  relative 
performance  of  the  various  configurations  rather  than  on  the 
absolute  performance  data  generated. 

The  study  indicates  that  performance  improvement  of 
about  15%  is  attained  either  by  reducing  the  preload  charge 
level  by  30%  or  by  reducing  the  MESG  errors  to  the  level 
specified  for  calibration  accuracy  (thus  not  including  er- 
rors due  to  repeatability  of  calibration  coefficients  across 
a shutdown).  Although  the  simulation  results  indicate  that 
better  navigation  performance  can  be  attained  by  reducing 
the  preload  charge  level,  this  change  was  not  incorporated 
into  the  MICRON  system  since  the  maximum  suspension  force 
that  can  be  generated  by  the  MESG  would  be  too  small.  The 
suspension  force  is  related  to  the  square  of  the  charge 
level;  consequently,  reducing  the  preload  charge  by  30% 
reduces  the  maximum  suspension  force  by  approximately  50%. 


Reducing  the  total  residual  MESG  errors  to  the  levels  speci- 
fied for  calibration  accuracy  is  an  available  design  option 
achievable  by  improving  total  instrument  performance  and/or 
improving  the  MESG  error  compensation  parameter  stability 
across  a power  cycle.  These  studies  present  the  performance 
improvement  which  can  be  expected  from  this  change.  The  ex- 
pected performance  gain  must  be  traded  against  the  other 
MICRON  targets  (cost,  reliability,  etc.)  prior  to  reaching 
a design  decision. 

2.3  VIBRATION  STUDIES 

The  MCSP  was  also  used  to  investigate  the  sensitiv- 
ity of  navigation  accuracy  to  correlated  linear  and  angular 
vibrations.  As  shown  in  Fig.  2.3-1,  the  vibro-pendulous 
motion  that  was  used  for  this  study  could  be  generated  by 
mounting  the  MICRON  system  on  the  end  of  a pendulum.  The 
three  parameters  which  characterize  this  motion  are:  LMA, 

the  length  of  the  moment  arm;  Bmax,  the  maximum  angular 
displacement,  and  u>p,  the  frequency  of  the  vibration.  A 
detailed  description  of  the  vibro-pendulous  motion  trajec- 
tory is  presented  in  Ref.  3. 

(•Mil 


Figure  2.3-1  Vibro-Pendulous  Motion  Model 

for  MCSP  Vibration  Studies 
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It  can  be  shown  that  the  maximum  acceleration  normal 
to  the  moment  arm  on  a body  under  the  vibro-pendulous  motion 
described  in  Fig.  2.3-1  is 


LMA  w 6 

p max 


(2.3-1) 


For  these  studies,  the  maximum  angular  displacement  was  de- 
fined in  terms  of  the  vibration  frequency  as 


max  2 

WP 


(2.3-2) 


where  k * 1 sec*2.  Because  of  the  definition  in  Eq.  (2.3-2), 
the  maximum  normal  acceleration  can  be  varied  independently 
of  the  frequency  of  vibration  by  changing  the  value  of  LMA. 

The  vibration  environment  studies  were  principally 
concerned  with  the  sensitivity  of  the  navigation  accuracy  to 
input  frequencies  close  to  the  fast  and  slow  loop  frequencies 
of  the  MICRON  mechanization.  The  fast  loop  operations  which 
consist  of  sampling  the  inertial  sensors  occur  at  64  times 
per  second  while  the  slow  loop  solution  of  the  navigation 
equations  is  completed  16  times  per  second.  Table  2.3-1 
presents  the  vibration  frequencies  and  maximum  normal  ac- 
celeration levels  that  were  used  in  the  MCSP  vibration 
studies.  Also  presented  in  the  table  is  the  value  of  the 
truth  loop  frequency  (TLF)  used  for  each  run.  As  described 
in  Ref.  3,  TLF  is  an  MCSP  input  parameter  that  controls  the 
frequency  of  solution  of  the  truth  model  equations  gener- 
ating the  reference  values  against  which  the  MICRON  naviga- 
tion outputs  are  compared. 

Figures  2.3-2  and  2.3-3  present  single-trial  navi- 
gation error  time  histories  for  30-second  MCSP  runs  for 
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TABLE  2.3-1 


VIBRATION  FREQUENCIES  AND  NORMAL 
ACCELERATIONS  FOR  MCSP  VIBRATION  STUDIES 


CASE 

VIBRATION 
FREQUENCY  (Hz) 

MAX  an 
(g’s) 

TLF 

Baseline 

0 

0 

64 

1 

120 

0.10 

256 

2 

120 

0.03 

256 

3 

90 

0.10 

256 

4 

90 

0.03 

256 

5 

GO 

0.10 

128 

6 

60  * 

0.03 

128 

7 

20 

3.00 

64 

8 

5 

64 

the  baseline  and  Case  1,  respectively.  For  all  the  vibration 
studies,  the  nominal  residual  error  parameters  were  used  to 
describe  the  MESG's.  As  illustrated  by  Fig.  2.3-3,  the  vibro- 
pendulous  environment  causes  oscillatory  velocity  errors. 
Because  of  the  oscillatory  nature  of  the  velocity  errors  in 
the  vibration  environment.  Table  2.3-2  (which  summarizes  the 
results  of  the  MCSP  vibration  studies)  presents  the  maximum 
errors  during  the  30  seconds  of  vibration  rather  than  the 
terminal  navigation  errors. 

The  results  of  the  MCSP  vibration  studies  indicate 
that  vibro-pendulous  input  motion  does  not  cause  any  unac- 
ceptable navigation  errors.  The  errors  caused  by  the  low 
amplitude,  high  frequency  vibrations  are  quite  small  rela- 
tive to  MICRON'S  performance  specifications  of  1 nm/hr 
position  error  growth  rate  and  5 ft/sec  velocity  error. 

The  lower  frequency,  high  g-level  vibrations  generated  siz- 
able navigation  errors.  Reference  10  presents  a plot  of  the 
vibration  amplitude  as  a function  of  the  frequency  to  which 
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Baseline  Single-Trial  Navigation  Errors 
for  MCSP  Vibration  Studies 


Figure  2.3-2 


Case  1 Single-Trial  Navigation  Errors 
for  MCSP  Vibration  Studies 


Not#:  Velocity  error#  ar#  owlllatlng  at  the  frequency  of  the  vlbro-pendulou#  motion. 


Figure  2.3-3 
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TABLE  2.3-2 


SINGLE-TRIAL  NAVIGATION  ERRORS  FOR 
MCSP  VIBRATION  STUDIES 


CASE 

MAX.  NORTH 
POSITION  ERROR 
(ft) 

MAX.  EAST 
POSITION  ERROR 
(ft) 

MAX.  NORTH 
VELOCITY  ERROR 
(ft/sec) 

MAX.  EAST 
VELOCITY  ERROR 
(ft/sec) 

Baseline 

-0.038 

0.460 

0.012 

0.054 

1 

-0.100 

0.423 

0.008 

0.036 

2 

-0.060 

0.368 

0.008 

0.032 

3 

-0.463 

0.790 

-0.039 

0.068 

4 

-0.161 

0.491 

-0.012 

0.042 

5 

-0.171 

0.545 

-0.033 

0.078 

6 

-0.062 

0.435 

0.014 

0.049 

7 

-78.193 

78.008 

-5.249 

5.875 

8 

-95.637 

95.514 

-10.000 

11.500 

the  MICRON  equipment  inside  its  shock  mounts  will  be  exposed 
The  maximum  amplitude  at  5 Hz  is  0.2g  so  the  MICRON  system 
will  see  very  little  of  the  motion  that  was  used  in  Case  8 
of  the  MCSP  vibration  studies.  The  resonant  frequency  of 
the  proposed  shock  mounts  is  20  Hz  (Ref.  10);  consequently, 
the  sensitivity  of  navigation  performance  to  this  frequency 
vibration,  as  indicated  by  Case  7 of  the  MCSP  vibration 
studies,  requires  additional  investigation. 


2.4  HIGH  ROLL  RATE  STUDIES 


The  final  MCSP  investigation  considered  the  effects 
on  MICRON  navigation  accuracy  of  severe  aircraft  maneuvers. 
These  maneuvers  included  "high-g"  linear  accelerations  and 
pullups  and  the  maximum  allowable  roll  rate  and  roll  ac- 
celeration as  delineated  by  the  MICRON  specification.  To 
isolate  the  errors  induced  by  the  high  roll  environment, 
two  trajectories  were  used  in  these  studies.  The  two  tra- 
jectories were  identical  except  that  the  roll  motion  envi- 
ronment was  included  in  only  one  of  the  trajectories. 


J | 


The  basic  trajectory  selected  was  Trajectory  4 de- 
scribed in  Ref.  3.  Figure  2.4-1  presents  an  altitude  profile 
of  this  30-second  trajectory.  On  a constant  northeast  head- 
ing, the  aircraft  undergoes  a 4g  linear  acceleration  to  a 
velocity  of  660  ft/sec  followed  by  a 5g  pullup  to  45° , 
and  a 6 second  climb.  The  climb  is  concluded  with  a 5g  turn 
down  maneuver  to  level  unaccelerated  flight.  During  the 
climb  maneuver  of  Trajectory  4,  the  total  velocity  magnitude, 
VT,  is  held  constant.  The  horizontal  velocity  VH  and  the 
vertical  velocity  h are  computed  at  each  iteration  of  the 
MCSP  Trajectory  Generator  from  the  aircraft  pitch  angle,  6. 
The  equations  which  determine  Vjj  and  h are 


V,p  cos  0 


(2.4-1) 


VT  sin  0 


(2.4-2) 


Table  2.4-1  presents  the  parameter  values  which  define  Tra- 
jectory 4 (consult  Ref.  3 for  a complete  description  of  all 
the  equations  used  by  the  MCSP  to  compute  the  reference  tra- 
jectory). 
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Figure  2.4-1 


Trajectory  4 - Acceleration 
and  Pullup 
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TABLE  2.4-la 

TRAJECTORY  4 INITIAL  VALUES 


Longitude 

X 

-72  degrees 

Vertical  Velocity 

h 

0 ft/sec 

Latitude 

L 

30  degrees  * 

Heading 

4- 

w/4  radians 

Altitude 

h 

0 feet 

Pitch  Angle 

e 

0 radians 

Total  Velocity 

VT 

0 ft/sec 

Roll  Angle 

♦ 

0 radians 

TABLE  2. 4 -lb 

TRAJECTORY  4 PARAMETER  VALUES  VERSUS  TIME 


Tim  (see) 

0-2 

2-7.5 

7.5-11 

11-15 

15-21 

21-2.1 

25-30 

VT 

(ft/sec) 

0 

% ♦/♦»<« 

660 

M0 

600 

660 

660 

*T  , 
(ft/ssc3) 

0 

120 

0 

0 

0 

0 

0 

♦ 

(radians) 

w/4 

»/4 

*/4 

t/4 

•/« 

»/4 

»/4 

♦ 

(radlans/seo) 

0 

0 

0 

0 

0 

0 

0 

0 

(radians) 

0 

0 

0 

•o  ♦ /*  « 

«/« 

•o  ♦ / • « 

0 

i 

(radlaaa/see) 

0 

i 

0 

0 

»/!• 

0 

-w/16 

0 

The  high  roll  rate  trajectory  (which  will  be  called 
Trajectory  4a)  superimposes  roll  motion  of  the  aircraft  on 
the  maneuvers  that  comprise  Trajectory  4.  The  roll  motions 
consist  of  a 1.2  second  roll  through  an  angle  of  140°  near 
the  completion  of  the  pull-up  maneuver  shown  in  Fig.  2.4-1 
and  a 1.2  second  roll  of  -140°  at  the  beginning  of  the  level- 
ing maneuver  which  concludes  the  climb  phase  of  the  trajec- 
tory. Table  2.4-2  presents  the  time  histories  of  the 
parameters  which  define  Trajectory  4a.  The  initial  condi- 
tions for  Trajectory  4a  are  identical  to  those  in  Table 
2.4-la  for  Trajectory  4. 


a 

0 
0 
Li 

01 

f 


L! 

0 


0 

0 
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TRAJECTORY  4a  - ROLL  RATE  ENVIRONMENT 
PARAMETER  VALUES  VERSUS  TIME 


Three  MCSP  runs  were  made  using  Trajectories  4 and 
4a  to  complete  the  high  roll  environment  studies.  Both 
trajectories  were  used  for  cases  employing  perfect  MESG's 
(no  residual  drift  rate  or  ARO  error)  to  determine  the  ef- 
fects of  mechanization  errors  induced  by  the  severe  environ- 
ments. For  the  high  roll  rate  trajectory,  a single  MCSP 
run  was  also  made  using  nominal  drift  rate  and  ARO  error 
(|e|  * 0.01  deg/hr,  |6y|  * 0.1  mrad)  residuals  to  describe 
the  MESG.  Figures  2.4-2  and  2.4-3  present  time  histories 
of  the  navigation  errors  during  Trajectory  4a  for  zero  and 
nominal  residual  errors  of  the  MESG,  respectively.  Table 
2.4-3  presents  the  terminal  navigation  errors  for  all  three 
cases . 

The  simulation  results  indicate  that  the  navigation 
errors  due  to  system  mechanization  (i.e.,  the  zero  MESG  er- 
ror runs)  are  acceptably  small  for  both  the  high  acceler- 
ation and  the  high  roll  rate  trajectories.  Thus,  the  fre- 
quencies chosen  for  sampling  the  inertial  instruments  and 
solving  the  navigation  equations  appear  to  be  adequate  for 
the  specified  angular  rate  environment.  The  differences 
between  the  Trajectories  4 and  4a  results  for  the  zero  MESG 
error  runs  indicate  that  the  solution— rate— induced  errors 
are  dependent  on  the  motion  to  which  the  navigation  system 
is  exposed. 

Introduction  of  nominal  residual  MESG  errors  into 
the  MCSP  yields  significantly  larger  position  and  velocity 
errors  for  the  high  angular  rate  trajectory  than  were  ob- 
served when  there  were  no  MESG  errors.  Because  of  the 
short  duration  of  this  trajectory  and  the  assumed  zero- 
valued initial  attitude  errors,  it  is  likely  that  the  navi- 
gation errors  are  due  principally  to  ARO  error.  (There 
has  not  been  sufficient  time  for  residual  drift  rate  to 
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Navigation  Errors  for  High  Roll  Rate 
Trajectory  (No  MESG  Errors) 


Figure  2.4-2 


Figure  2.4-3  Single-Trial  Navigation  Errors  for 

High  Roll  Rate  Trajectory  (Nominal 
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TABLE  2.4-3 


TERMINAL  NAVIGATION  ERRORS  FOR 
MCSP  HIGH  ROLL  RATE  STUDIES 


CASE 

THSUNAL  NORTH 
POSITION  DWOR 

(ft) 

TDWINAL  EAST 
POSITION  DOOR 

(ft) 

TOW  INAL  NORTH 

va/jcm  0WOR 

(ft/BOC) 

TOW  INAL  EAST 
vnrcm  men 

(ft/sec) 

Trajectory  4 
(No  MESO  Error*) 

•0*226 

-0.265 

-0.164  x 10~3 

-3.47  x 10": 

Trajectory  4a 
(No  MESG  Errors) 

-0.262 

-0.051 

-1.93  * 10"3 

-3.80  x 10": 

Trajectory  4a 
(Nominal  MESG  Errors) 

1.367 

0.193 

106.0  x 10"3 

33.0  x 10"; 

propagate  into  navigation  errors.)  ARO  errors  during  the 
aircraft  maneuvers  cause  a misresolut ion  of  the  sensed  accel- 
eration resulting  in  a velocity  error.  The  single-trial  MCSP 
results  for  the  high  roll  rate  trajectory  and  nominal  MESG 
errors  are  consistent  with  MICRON  performance  targets  (Ref.  5) 
The  conclusion  to  be  drawn  from  these  studies  is  that  high 
roll  rates  and  accelerations  do  not  prevent  MICRON  from 
achieving  its  specified  performance  goals. 
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COVARIANCE  ERROR  ANALYSIS  STUDIES 


The  MICRON  Computer  Simulation  Program  (MCSP)  which 
was  exercised  to  generate  the  results  presented  in  Chapter  2 
is  a direct  simulation  of  the  system  navigation  equations. 

It  provides  the  navigation  accuracy  corresponding  to  the 
deterministic  set  of  conditions  defined  by  the  program  in- 
put parameters.  To  conduct  statistical  performance  analy- 
sis studies  using  the  MCSP,  it  would  be  necessary  to  make 
a large  number  of  computer  runs  using  input  error  sources 
provided  by  random  number  generators  with  appropriate  sta- 
tistics. Computing  the  ensemble  statistics  of  the  results 
of  these  runs  would  yield  the  desired  statistical  perform- 
ance prediction.  (This  technique  is  called  Monte  Carlo 
analysis. ) 


"Covariance  Analysis"  is  a well-known  alternative 
to  Monte  Carlo  techniques  (Ref.  11).  Using  covariance 
analysis,  a single  computer  run  is  used  to  generate  a 
prediction  of  the  MICRON  system  performance  over  an  en- 
semble of  similar  missions.  Covariance  analysis  evalua- 
tion of  inertial  navigation  system  performance  is  a com- 
monly accepted  technique  which  has  been  used  successfully 
during  all  stages  of  navigation  system  development  pro- 
grams (Refs.  12  and  13).  By  reviewing  the  mechanization 
equations  and  developing  models  for  MESG  residual  errors, 
this  chapter  will  describe  the  preliminary  covariance 
error  analysis  program  generated  by  TASC  for  use  in  MICRON 
analysis  studies.  Since  at  this  stage  in  MESG  instrument 
development,  there  has  not  been  sufficient  data  gathered 
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to  determine  valid  statistical  models  for  the  residual  drift 
rate  and  ARO  errors,  it  was  assumed  for  the  preliminary  co- 
variance  analysis  that  the  errors  were  due  to  miscalibra- 
tion  of  some  of  the  simpler  drift  rate  and  ARO  model  terms. 
The  mission  assumed  for  these  early  analyses  was  navigation 
in  a stationary  laboratory  environment.  Comparison  of  the 
covariance  results  with  lab  test  data  will  provide  insights 
allowing  improvement  of  the  MESG  residual  error  models  dur- 
ing future  studies. 

During  Phase  2 of  MICRON  development,  after  the  MESG 
design  has  been  fixed,  it  is  suggested  that  a Sensor  Data 
Bank  be  established  to  collect  data  on  a large  number  of 
MESG's.  This  data  should  then  be  used  to  develop  statisti- 
cal models  that  accurately  describe  the  residual  drift  rate 
and  ARO  errors  of  the  MESG.  These  improved  models  can  be 
used  to  refine  the  error  analysis  program,  which  was  devel- 
oped as  part  of  the  current  effort.  Expanding  the  mission 
scenario  capability  will  allow  the  evaluation  of  MICRON 
performance  in  any  proposed  application  of  interest. 

The  model  of  MICRON  system  performance  to  be  gen- 
erated during  these  covariance  analysis  studies  is  the  kind 
of  system  description  required  for  using  modern  data  proces- 
sing techniques  (i.e.,  Kalman  filtering)  to  integrate  MICRON 
with  those  navaids  that  are  available  in  a modern  avionics 
suite.  Thus,  in  addition  to  providing  system  performance 
predictions,  development  of  covariance  error  analysis  tools 
for  MICRON  will  facilitate  future  development  of  aided  in- 
ertial system  mechanizations. 

Section  3.1  describes  the  mechanization  proposed  for 
use  with  the  MICRON  system.  The  error  equations  correspond- 
ing to  this  mechanization  are  derived  in  Section  3.2,  while 
the  error  analysis  computer  program  and  some  typical  covar- 
iance results  are  presented  in  Section  3.3. 


3.1  OVERVIEW  OF  MECHANIZATION 


Prior  to  describing  the  error  equations  for  MICRON, 
it  is  useful  to  discuss  briefly  the  navigation  mechanization 
equations  which  have  been  implemented.  Although  M I CRON  is 
a strapdown  system,  the  MESG  rotors  remain  fixed  in 
space  allowing  determination  of  vehicle  attitu  e r 
from  the  gyros.  Consequently,  the  typical  strap  ^ 

anization  problems  associated  with  Integra  ng 
do  not  exist  in  MICRON.  The  fact  that  the  MICRON  acceler- 
ometers are  mounted  directly  to  the  airframe  Presents  a 
classical  strapdown  problem  of  properly  integrating  the 
accelerometer  outputs  to  avoid  the  corrupting  effects 
vehicle  attitude  rates. 


The  MICRON  mechanization  (Ref.  14)  uses  three 
mate  frames:  the  MICRON  or  body  frame  (designated  the , M- 

frame)  which  is  fixed  to  the  vehicle,  the  S-frame  (spin)  de 
fined  by  the  rotor  spin  axes  of  the  two  MESG’s,  and  a local- 
level,  free  azimuth  navigation  (N)  frame.  The  M-frame  ha 
the  input  axis  of  one  of  the  accelerometers  along  each  o 
its  three  axes.  Also,  the  case-fixed  coordinate  frames^ 
both  the  MESG's  are  aligned  with  the  M-  rme. 

MESGA  multisensor  replaces  the  accelerometers  in  the  f 

.h  M frame  will  likely  be  defined  by  the 
MICRON  system,  the  M-frame  w 

case-fixed  axes  of  the  MESGA.) 


The  defining  axes  of  the  S-frame  are  Ij , X2'  a"d . 

Y x - where  1(  is  the  direction  of  the  rotor  spin  axis 
*1  ‘ of 2 the  i«>  JeSG.  This  frame,  which  is  inertially  fixed 

xcept  or  the  drift  o,  the  RSA's.  is  not  orthogonal  because, 
although  the  gyros  are  spun  up  approximately  90  apart,  the 

two  RSA's  are  not  torqued  and  consequently  -constra 
to  be  orthogonal.  The  fact  that  the  S-frame  is  not  orthogo 
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does  not  present  any  severe  difficulties  as  long  as  the  ap- 
propriate transformations  account  for  the  nonorthogonality. 

Figure  3.1-1  presents  a flowchart  describing  the 
MICRON  mechanization  equations.  There  are  two  solution  loops 
associated  with  this  implementation:  an  inner  loop  or  "fast 

cycle"  which  performs  time-critical  computations  64  times  a 
second  and  an  outer  loop  or  "slow  cycle"  which  executes  the 
remainder  of  the  navigation  calculations  at  16  times  per 
second  (Ref.  14).  The  strapdown  accelerometer  problem  in 
MICRON  is  solved  by  the  fast  cycle  operations.  During  each 
fast  cycle,  the  accelerometer  outputs  are  projected  onto  the 
defining  axes  of  the  S-frame  (determined  from  the  MESG  atti- 
tude readout),  transforming  the  sensed  acceleration  during 
the  fast  cycle  interval  into  the  S-frame.  This  process  cre- 
ates a "mathematical"  accelerometer  triad  (one  accelerometer 
along  each  axis  of  the  S-frame)  that  is  essentially  iner- 
tially  fixed  and  decoupled  from  vehicle  attitude  rates  (as 
would  be  the  case  on  a space  stable  gimballed  platform). 

The  integrated  accelerations  in  the  S-frame  are  then  used 
in  the  slow  loop  for  the  remainder  of  the  navigation  compu- 
tations . 

The  integrated  acceleration  in  the  S-frame  is  in  er- 
ror due  to  the  effects  of  accelerometer  errors  and  MESG  atti- 
tude readout  (ARO)  errors.  Compensation  for  these  effects 
using  the  accelerometer  and  ARO  error  model  coefficients  is 
performed  during  the  slow  cycle.  Following  the  compensa- 
tion of  the  integrated  acceleration  in  the  S-frame,  the 
velocity  change  is  transformed  to  the  local-level  N-frame 
where  position  and  velocity  updates  occur  as  shown  in  Fig. 

3.1-1. 

The  transformation  from  the  S-frame  to  the  N-frame, 
Cs,  is  computed  as  a function  of  the  RSA  directions  coordin- 
ated in  the  N-frame.  (Since  the  axes  of  the  S-frame  are 
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* 
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not  orthogonal,  is  not  the  familiar  direction  cosine 

matrix  that  relates  orthogonal  coordinate  systems.)  It  is 

not  possible  to  measure  the  RSA  directions  in  the  N-frame 

since  the  MESG  attitude  readout  is  relative  to  the  gyro  case 

which  is  mounted  in  the  M-frame.  Consequently,  estimates 

of  the  RSA  vectors  in  the  N-frame  are  propagated  during  the 

slow  loop  using  the  MESG  drift  rate  model  and  the  angular 

rate  of  the  N-frame  relative  to  inertial  space  due  to  earth 

rate  and  vehicle  motion.  The  initial  values  of  the  N-frame 

estimates  of  the  RSA  vectors  (which  define  the  initial  value 

of  C«)  are  determined  by  the  MICRON  alignment  procedure  (see 
N 

Ref.  14). 

The  one  measurement  that  can  be  used  to  update  the 
estimates  of  the  N-frame  RSA  direction  vectors  is  the  angle 
between  and  defined  by  ^ ^ (the  dot  product  is  a 

scalar  and  thus  independent  of  coordinate  system  in  which  it 
is  computed).  In  the  MICRON  mechanization,  the  redundant 
axis  control  updates  the  estimate  of  the  N-frame  RSA  vector 
for  the  secondary  MESG  (that  gyro  whose  RSA  is  initially 
horizontal)  so  that  the  angle  between  the  N-frame  RSA  vec- 
tor estimate  tracks  the  measured  angle  in  a feedback  config- 
uration. Redundant  axis  control,  as  will  be  discussed  in 
the  error  analysis  development,  forces  the  error  introduced 
by  resolving  the  sensed  acceleration  vector  from  the  M-frame 
(where  it  is  measured)  to  the  N-frame  (where  the  navigation 
equations  are  solved)  through  the  nonorthogonal  S-frame  to 
be  the  small  vector  angle  rotation  that  is  typical  of  all 
inertial  systems.  (Reference  9 presents  a discussion  of 
small  vector  angle  rotations.)  The  next  section  will  de- 
velop the  error  equations  that  describe  the  MICRON  system. 
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3.2  ERROR  EQUATION  DEVELOPMENT 

The  dynamic  propagation  of  basic  error  sources  into 
INS  system  errors  is  well  known  (Refs.  9 and  15).  It  has 
been  shown  that,  for  any  dynamically  exact  implementation 
of  the  terrestrial  navigation  equations,  the  propagation 
of  inertial  instrument  errors  and  gravity  model  errors  can 
be  viewed  independently  of  the  particular  mechanization  that 
actually  appears  in  the  system  (Ref.  16).  Consequently,  the 
error  equations  can  be  written  in  any  convenient  coordinate 
system  as  long  as  care  is  taken  to  properly  identify  the  spe- 
cific driving  error  sources  for  the  system  and  mechanization 
being  studied. 

In  addition  to  instrument  induced  dynamic  errors  which 
are  independent  of  the  particular  mechanization,  there  also 
exist  implementation  errors  (caused  by  computer  roundoff,  num- 
erical inaccuracy,  and  solution  rates)  which  will  depend  on 
the  particular  mechanization  employed.  The  covariance  error 
analysis  technique  is  not  well  suited  for  investigating  the 
effects  of  implementation  errors;  and  thus,  they  will  not  be 
considered.  The  MCSP,  which  duplicates  the  MICRON  mechaniza- 
tion, has  been  used  to  investigate  these  errors  (Ref.  3). 

The  philosophical  approach  to  error  equation  devel- 
opment taken  in  Refs.  9 and  15  requires  the  Identifica- 
tion of  three  coordinate  frames,  ideally  coincident  but 
actually  differing  by  small  angles  due  to  navigation  errors. 
These  frames  are:  the  inertial  (often  called  the  platform) 

frame  which  is  an  inertially  stabilized  frame  where  the 
accelerometer  outputs  are  measured,  the  computation  frame 
in  which  the  navigation  equations  are  solved,  and  the  true 
frame  in  which  it  is  desired  to  mechanize  the  navigation 
equations.  Summarizing  the  results  of  the  derivation,  the 
dynamic  equations  of  INS  error  propagation  can  be  written  as 
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£ («£L)  ' - “el  x siL  + 6^L  (3'2"1) 

£ («vL)  - - (a£lj  ♦ s>el)  x s^L  - iL  x 


_d_ 

dt 


/ L L\  ,L  . L 

= - ^IE  + wELJ  x i + £ 


(3.2-3) 


where 

6r  * INS  position  error 

6v  = INS  velocity  error 

£ = the  small  vector  angle  rotation  such  that 

if  the  computation  frame  were  rotated  about 
the  £ direction  by  an  angle  \p  it  would  be 
brought  into  coincidence  with  the  inertial 
(or  platform)  frame. 

w = angular  rate  of  the  L-frame  with  respect 
EL  to  (w.r.t.)  an  earth-fixed  axis  set 

<*)..„  = angular  rate  of  an  earth-fixed  axis  set 
w.r.t.  inertial  space 

6g  * error  in  calculation  of  gravity 
= error  in  acceleration  sensing 
e = gyro  drift  rate 


The  superscript  "LM  refers  to  the  coordinate  frame  in  which 
the  error  analysis  is  to  be  performed.  The  driving  error 
terms  6g,  6f  and  £ depend  on  the  particular  system  and  mech- 
anization being  studied  and  will  be  defined  for  the  MICRON 
system  in  the  remainder  of  this  section. 


The  MICRON  error  equations  will  be  written  in  a north- 
slaved  local-level  coordinate  frame.  Although  MICRON  is  mech- 
anized using  a local-level  wander  azimuth  system  (since  it 
does  not  exhibit  the  polar  singularity  observed  in  north- 
slaved  systems),  it  is  desirable  to  write  the  error  equations 
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in  a north-slaved  frame  since  several  error  sources  (gravity 
errors  and  initial  condition  errors)  are  most  easily  modeled 
in  a north-slaved  frame.  Also,  output  errors  (navigation 
position  and  velocity)  are  usually  cited  in  this  frame. 

Thus,  for  the  MICRON  error  equation  development,  the  L- 
frame  used  in  Eqs.  (3.2-1)  to  (3.2-3)  will  be  the  local- 
level  north  slaved  (NED)  coordinate  system  with  its  x-axis 
north,  y-axis  east,  and  z-axis  down.  Section  3.2.1  will 
present  the  form  of  the  velocity  error  driving  terms  [«£ 
and  6f  in  Eq.  (3.2-2)]  for  the  MICRON  system  while  Section 
3.2.2  will  derive  the  MICRON  gyro  drift  rate  expression 
which  drives  the  equation  [Eq.  (3.2-3)]. 


3.2.1  MICRON  Gravity  Calculation  and 
Acceleration  Sensing  Errors 


Accelerometers  measure  specific  force,  f.,  defined  as 


f » a - £ 


(3.2-4) 


where  a is  the  inertial  acceleration  and  £ is  the  gravity  vec- 
tor. Since  the  inertial  navigation  computations  of  position 
and  velocicy  (which  consist  of  integrating  Newton's  second 
law  of  motion)  require  knowledge  of  the  inertial  accelera- 
tion, a,  it  is  necessary  that  the  system  mechanization  esti- 
mate £ using  a gravity  model  and  the  indicated  vehicle 
position.  While  avoiding  a detailed  discussion  of  the  grav- 
ity model  that  is  implemented  in  MICRON,  it  is  possible  to 
identify  two  classes  of  gravity  errors,  «£,  which  drive  the 
error  equations.  These  classes  of  errors  are:  (1)  errors 

due  to  evaluating  the  gravity  model  at  the  erroneous  navi- 
gator-indicated position,  and  (2)  errors  caused  by  a dif- 
ference between  the  true  gravity  and  the  gravity  model  eval- 
uated when  position  is  known  exactly. 
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It  can  be  shown  that  in  the  L-frame,  the  gravity 
error  due  to  the  INS  position  error  can  be  expressed  as 


& * 


(3.2-5) 


where  g is  the  nominal  value  of  gravity  and  R is  the  radius 
of  the  earth.  Inserting  Eq.  (3.2-5)  into  Eq.  (3.2-2)  and 
combining  with  Eq.  (3.2-1)  demonstrates  the  classical  oscil- 
lation of  horizontal  position  error  at  the  Schuler  frequency 
w = /So/R  and  the  instability  of  vertical  channel  navi- 
gator errors.  The  vertical  channel  is  stabilized  using  an 
altimeter  and  feeding  the  difference  between  the  INS  vertical 
position  and  the  altimeter  output  back  through  appropriate 
gains  to  the  inputs  of  the  vertical  position  and  velocity 
integrators.  The  altimeter  output  is  not  used  in  the  gravity 
calculation  for  MICRON.  This  stabilization  technique  intro- 
duces altimeter  errors  into  the  error  equations  for  the 
vertical  channel.  Appendix  A presents  a derivation  of  the 
vertical  channel  error  equations  resulting  from  altimeter 
damping. 

Gravity  errors  due  to  differences  between  the  true 
gravity  and  the  modeled  gravity  when  the  position  is  known 
exactly  are  called  gravity  anomalies  and  deflections  of  the 
vertical.  Gravity  anomalies  are  errors  in  the  magnitude 
of  the  model  estimate  while  deflections  of  the  vertical  are 
errors  in  the  direction  of  the  gravity  vector.  TASC  has 
conducted  extensive  studies  directed  both  at  developing 
statistical  descriptions  of  gravity  model  errors  (Ref.  17) 
and  investigating  the  effects  of  these  errors  on  INS  per- 
formance accuracy  (Ref.  18).  For  a moderate  accuracy  iner- 
tial system  such  as  MICRON,  gravity  anomalies  and  deflections 
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of  the  vertical  are  not  expected  to  be  dominant  error  sources. 
Consequently,  these  fig  error  terms  were  not  included  in  the 
preliminary  MICRON  error  equations. 


The  derivation  of  the  acceleration  sensing  errors, 


fifL,  which  drive  the  velocity  error  equation  [Eq.  (3.2-2)] 


requires  careful  review  of  the  MICRON  mechanization  to  in- 
sure that  the  errors  are  properly  transformed  from  the 
coordinate  frame  in  which  they  occur  into  the  L-frame  where 
the  error  equations  are  being  written.  It  will  be  shown 
that  the  MESG  Attitude  Readout  errors  drive  the  error  equa- 
tions in  a similar  manner  to  accelerometer  misalignments. 


To  avoid  confusion  caused  by  the  large  number  of  dif- 
ferent coordinate  frames,  it  is  useful  to  associate  the  co- 
ordinate frames  in  the  MICRON  mechanization  as  described  in 
Section  3.1  with  the  three  coordinate  frames  that  formed  the 
basis  of  the  error  equation  development.  The  MICRON  S-frame 
defined  by  the  spin  axis  vectors  of  the  two  MESG's  corresponds 
to  the  inertial  (platform)  frame  while  the  local -level,  azi- 
muth wander  N-frame  can  be  viewed  as  the  computation  frame. 

The  true  frame  is  that  frame  in  which  it  is  desired  to  mech- 
anize the  MICRON  navigation  equations  and  thus  is  defined 
by  the  actual  vehicle  position.  Although  it  does  not  cor- 
respond to  one  of  the  frames  used  to  derive  the  error  equa- 
tions, the  body  fixed  M-frame  (see  Section  3.1)  is  the 
coordinate  system  in  which  all  the  inertial  sensors  are 
mounted  and  thus  is  the  frame  where  the  instrument  errors 
actually  occur. 


As  discussed  in  Section  3.1,  the  MICRON  mechanization 


solves  the  strapdown  accelerometer  problem  by  projecting  the 
accelerometer  outputs,  fM,  onto  ^ , £2  and  yx  x y2  creating 
a mathematical  accelerometer  triad  in  the  inertially  stable 
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S-frame.  The  specific  force  in  this  S-frame,  fS,  can  be 
related  to  the  accelerometer  outputs  (i.e.,  specific  force 
in  the  M-frame),  , by 


(3.2-6) 


where 


I’m 

' M 

( M 

m\1t 

Ui 

: *2 

Ui 

x l2)\ 

(3.2-7) 

The  acceleration  sensing  errors  relative  to  the  platform  frame 
(which,  for  MICRON,  is  the  S-frame)  determine  the  system  driv- 
ing errors  6f  in  Eq.  (3.2-2). 


Perturbing  Eq.  (3.2-6)  yields 
6fS  = C®  (6fM  + 6C®  fM) 


6fM  is  the  accelerometer  instrument  error  and  6C^  is  the  er- 
ror  introduced  by  the  MESG  ARO  errors  in  constructing  the 
mathematical  accelerometer  triad  along  the  axes  of  the  S- 
frame.  The  errors  associated  with  the  accelerometer  are 
scale  factor  and  misalignments,  random  bias,  and  a first 
order  Markov  noise  error.  Including  these  terms,  it  is  pos- 
sible to  express  the  accelerometer  instrument  errors,  6f^,  as 

6fM  = K f M + 7b  + 7r  (3.2^9) 

The  matrix  K models  scale  factor  and  misalignment  errors 
whiie  7fe  and  7r  represent  the  bias  and  Markov  noise  compon- 
ents of  the  accelerometer  errors. 
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The  MESG  ARO  errors  that  remain  following  compensa- 
tion cause  the  readout  position  of  the  rotor  spin  axis  (RSA) 
to  be  misaligned  from  its  true  position.  Thus,  projecting 
the  accelerometer  data  onto  these  misaligned  RSA  vectors  will 
result  in  an  effective  misalignment  of  the  mathematical  ac- 
celerometer triad.  It  is  this  error  source  which  is  modeled 
by  the  6C®  fM  term  in  Eq.  (3.2-8).  By  perturbing  Eq.  (3.2-7) 
and  retaining  only  first  order  terms,  it  is  possible  to  com- 

g 

pute  60^  as 

{Cm  ' [5i“  : «x2M  : (x?  * si“)  + K * x2)] 

(3.2-10) 

where  Sy^  and  6^®  represent  the  residual  ARO  error  of  the 
primary  and  secondary  MESG’s,  respectively.  Using  Eq.  (3.2-10) 
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(3.2-11) 


The  form  of  the  statistical  models  describing  the 
residual  ARO  errors,  6^,  needed  to  allow  the  use  of  Eq. 
(3.2-11)  in  the  MICRON  error  analysis  are  not  yet  well  de- 
fined (due  to  insufficient  MESG  data).  For  the  preliminary 
development  of  the  error  analysis  computer  program  that  was 
performed  as  part  of  the  current  effort,  the  residual  ARO 
errors  were  assumed  to  be  due  to  miscalibration  of  some  of 
the  principal  terms  in  the  ARO  compensation  model  (Ref.  3). 
During  future  MICRON  analysis  studies,  when  the  MESG  design 
has  been  frozen  and  there  is  data  available  on  a large 
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number  of  gyros,  it  is  suggested  that  investigations  be  pur- 
formed  to  develop  valid  models  to  describe  the  residual 
MESG  instrument  errors. 


The  acceleration  sensing  errors  described  by  Eq.  (3.2-8) 
are  coordinat ized  in  the  S-frame.  Since  it  has  been  decided 


to  perform  the  MICRON  error  analysis  in  the  L-frame,  it  is 
necessary  to  transform  6f^  into  the  L-frame  before  it  is  used 


in  Eq.  (3.2-2).  Thus 


<4  c*  («“  * «$  f“) 

cjj  (fifM  + 6C§  fM) 


(3.2-12) 


The  transformation  Cjj,  which  relates  the  body- fixed  M-frame 
to  the  local-level,  north  slaved  L-frame,  will  be  a function 
of  the  vehicle  attitude  and  consequently  will  depend  on  the 
trajectory  for  which  the  error  analysis  is  performed. 


Summarizing  by  substituting  Eqs.  (3.2-9)  and  (3.2-11) 
into  Eq.  (3.2-12),  it  is  possible  to  express  the  accelera- 
tion sensing  errors  (accelerometer  and  ARO  errors)  that  drive 
the  velocity  error  equation  [Eq.  (3.2-2)]  for  the  MICRON  sys- 


tem as 


I r ♦ V.  ♦ V.  ♦ 
i — — o — r 


«ii“  ♦ 


I . « 


!**  *xi 


(xi  * iuf 


tf  * Xi)’ 


(3.2-13) 


M M 

where  the  models  for  the  residual  ARO  error  6^.  and  fijfg 


remain  to  be  defined.  The  models  that  were  used  during  the 
development  of  the  error  analysis  computer  program  will  be 
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presented  in  Section  3.3.1.  Appendix  B will  delineate  the 
system  state  vector  and  its  dynamics  matrix  that  described 
the  MICRON  navigation  errors. 


3.2.2  Gyro  Drift  Error  Propagation 

In  the  previous  section,  the  error  sources  which  drive 
the  velocity  error  equation  for  the  MICRON  system  were  identi- 
fied. To  complete  the  MICRON  error  analysis  equations,  it  is 
necessary  to  derive  the  driving  error  sources  for  the  £ angle 
©quation  [Eq.  (3.2— 3)3 . Ref.  9,  which  presents  a detailed 
derivation  of  the  ^ angle  equation,  views  the  vector  angle 
rotation,  as  descriptive  of  the  error  introduced  in  trans- 
forming a vector  (in  particular,  the  sensed  acceleration  vec- 
tor) from  the  stabilized  platform  frame  to  the  computation 
frame.  This  error  is  caused  by  an  error  in  the  computer's 
estimate  of  the  transformation  from  the  platform  frame  to 
the  computation  frame. 

For  the  MICRON  mechanization,  the  S-frame  defined  by 
the  MESG  RSA  direction  vectors  has  been  identified  as  the 
platform  frame  while  the  computation  frame  is  the  local-level 
wander  azimuth  N-frame  defined  by  the  computed  vehicle  posi- 
tion and  wander  angle.  Thus,  ^ describes  the  error  in  the 
estimate  of  S-frame  basis  vectors  (i.e.,  the  RSA  vectors) 
in  N-frame  coordinates  since  it  is  the  estimates  of  the  spin 
vectors  (in  the  N-frame)  which  define  Cg  . The  estimation 
process  assumes  that  the  RSA's  are  fixed  relative  to  iner- 
tial space,  except  for  the  redundant  axis  control  and  the 
modeled  gyro  drift.  Consequently,  residual  (unmodeled) 
drift  of  the  MESG  spin  axes  relative  to  inertial  space  is 
the  only  system  error  source  dynamically  driving 

The  redundant  axis  control  constrains  the  estimate 
of  the  spin  axis  vector  for  the  secondary  MESG  (the  primary 
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gyro's  RSA  is  initially  vertical  while  the  RSA  of  the  sec- 
ondary MESG  is  horizontal  and  forward  along  the  aircraft 


longitudinal  axis)  in  N-frame  coordinates,  SN2. 


As  is  dem- 


onstrated in  Appendix  C,  the  redundant  axis  control  adjusts 
SN2  in  the  SN1  - SN2  plane  so  that  the  scalar  product 

-N1  ' -N2  ec>uals  Xci  ' Xc2  which  is  determined  from  the 
compensated  readout  of  the  RSA  vectors.  Consequently,  for 
the  secondary  gyro  only  that  residual  drift  which  moves  Sn2 
out  of  the  SN1  - SN2  plane  will  affect  the  £ angle. 

SPj  and  6P2  are  now  defined  as  the  residual  drift 
rate  vectors,  coordinat ized  in  the  MESG  case  coordinates  (M- 
frame),  of  the  primary  and  secondary  gyros,  respectively. 

It  is  then  possible,  using  the  constraint  imposed  by  redund- 
ant axis  control,  to  determine  an  expression  for  e [Eq. 
(3.2-3)]  in  terms  of  6PX  and  6P2.  This  expression  is 


(3.2-14) 


M 

since  residual  drift  rate  of  the  primary  gyro  in  the  ^-di- 
rection does  not  lead  to  any  error;  and,  due  to  redundant 
axis  control , only  the  residual  drift  rate  about  the  ^-di- 
rection of  the  secondary  MESG  causes  angle  errors. 

Since  the  error  analysis  is  to  be  performed  in  the 
L- frame , it  is  necessary  to  transform  the  drift  rate  vector 
in  Eq.  (3.2-14)  into  the  L-frame.  Thus 

£L  - [(*  - II  if)  {£l  + I?  if  SE2]  (3.2-15) 

Summarizing  the  results  of  this  section,  Fig.  3.2-1 
presents  the  MICRON  error  propagation  equations  in  vector- 
matrix  form.  The  next  section  will  describe  the  computer 
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program  which  was  used  for  the  MICRON  error  analysis  studies 
and  will  present  some  of  the  preliminary  results  that  were 

generated . 


3.3  ERROR  ANALYSIS  COMPUTER  PROGRAM 


] 


To  allow  prediction  of  the  MICRON  system  performance, 
the  covariance  equations  corresponding  to  the  error  dynamics 
derived  in  Section  3.2  have  been  coded,  forming  the  error 
analysis  computer  program.  The  preliminary  version  of  this 
program  (which  is  described  in  this  section)  is  specialized 
to  a stationary  laboratory  environment  with  discrete  atti- 
tude changes  and  assumes  that  the  residual  drift  rate  and 
ARO  errors  are  due  to  miscalibration  of  terms  appearing  in 
the  compensation  models.  During  future  MICRON  error  covar- 
iance simulation  studies,  it  is  suggested  that  the  computer 
program  be  refined  to  incorporate  the  trajectories  of  planned 
MICRON  applications  and  to  include  more  accurate  descriptions 
of  the  residual  instrument  errors  as  they  become  available. 


The  state  vector  description  of  MICRON  system  errors 
used  in  the  error  analysis  computer  program  consists  of  45 
elements.  Nine  of  these  states  (3  position  errors,  3 vel- 
ocity errors,  and  3 ^-angles)  represent  the  system  accuracy 
while  the  remaining  36  states  model  the  error  sources  which 
drive  the  system  states.  Table  3.3-1  presents  a list  of  the 
error  sources  included  in  the  covariance  analysis  program 
and  the  number  of  states  which  are  required  to  model  these 


errors. 

As  was  presented  in  Section  3.2,  the  residual  drift 
rate  and  ARO  errors  of  the  individual  MESG’s  propagate  into 
system  errors  as  a function  of  the  RSA  directions.  Also, 
some  of  the  models  used  to  describe  residual  instrument 
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TABLE  3.3-1 
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errors  depend  on  the  RSA  direction.  Consequently,  the  dynam- 
ics of  many  of  the  system  driving  states  are  related  to  the 
motion  of  the  inertially-f ixed  RSA’s  which,  for  the  station- 

Dary  laboratory  environment,  move  at  earth  rate  relative  to 

the  L-frame  where  the  error  analysis  is  to  be  performed.  The 
fact  that  the  gravity  vector  is  constant  and  oriented  along 
one  of  the  axes  of  the  local-level  L-frame  has  been  used  to 
reduce  the  number  of  states  required  to  describe  the  MICRON 

D 
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MICRON  ERROR  ANALYSIS  PROGRAM  STATES 


DESCRIPTION 


NUMBER  OF 
INDIVIDUAL 
SOURCES 


NUMBER  OF 
ADDED  STATES 


MICRON  Position  Error 

MICRON  Velocity  Error 

tji-Angle  Misalignments 

Accelerometer  Bias  Errors 

Accelerometer  Random  Errors 

Accelerometer  Scale  Factor 
Errors 

Altimeter  Error 
ARO  Scale  Factor  Errors 
Case-Fixed  Bias  Drifts 
Random  Drift 

Acceleration  Sensitive  Drift 
Miscalibration 

Acceleration  Sensitive  Drift 
Miscompensation  Due  to 
Accelerometer  Error 

TOTAL 


3 

3 

3 


6 

(3  for  each  MESG) 
6 

(3  for  each  MESG) 


1 

1 

42 


3 

3 

3 

3 

3 

1 

1 

IS 

9 

3 

1 

0 

45 


error  sources.  When,  during  future  studies,  vehicle  motion 
is  added  to  the  error  analysis  program,  it  is  likely  that 
additional  states  will  be  required  to  properly  account  for 
the  effects  of  the  MICRON  error  sources. 

Section  3.3.1  describes  the  instrument  residual  er- 
ror models  that  have  been  used.  (Appendix  B presents  the  state 
vector  and  dynamics  matrix.)  Section  3.3.2  presents  the  com- 
puter program  description  while  Section  3.3.3  presents  some 
of  the  results  of  the  checkout  runs  that  were  performed  with 
the  program. 


3.3.1  Residual  Instrument  Errors 

A brief  description  of  the  sensor  errors  that  have 
been  included  in  the  error  analysis  program  is  presented. 

For  drift  rate  and  ARO  errors,  the  selected  models  have  been 
determined  by  ascribing  the  residual  errors  to  miscalibra- 
tion  of  dominant  compensation  terms.  During  future  MICRON 
Analysis  efforts,  it  is  suggested  that  test  data  from  a large 
number  of  MESG's  be  evaluated  to  determine  valid  statistical 
descriptions  of  residual  drift  rate  and  ARO  errors. 

Accelerometer  Errors  — Three  types  of  errors  have 
been  associated  with  the  accelerometers  for  the  MICRON  sys- 
tem model  used  in  the  error  analysis  program.  These  are 
bias,  scale  factor  and  random  errors.  The  random  acceler- 
ometer errors  are  modeled  by  a first  order  Markov  process. 
Although  the  error  equations  are  written  in  the  L-frame,  the 
accelerometers  are  physically  mounted  in  the  M-frame  (see 
Section  3.2).  Consequently,  the  program  input  parameters 
which  describe  the  accelerometer  errors  are  defined  rela- 
tive to  the  M-frame,  and  the  error  analysis  program  must 
compute  the  equivalent  errors  in  the  L-frame  for  use  in  the 
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4- 


MICRON  error  equations.  It  is  thus  possible  to  express  the 
total  accelerometer  error  vector  in  the  L-frame,  v£,  as 


+ V M 
— r 


+ 6C , 


(3.3-1) 


where  vj1  represents  the  bias  errors . 


models  the  random 


acclerometer  errors,  and  6C, 


is  a diagonal  matrix  of  scale 


factor  errors  for  the  three  accelerometers. 

Altimeter  Errors  — The  typical  errors  that  occur 
in  altimeters  are  bias  and  scale  factor.  For  the  stationary 
laboratory  environment,  the  altitude  is  constant  and  the  scale 
factor  error  in  the  altimeter  is  constant.  Thus,  it  is  pos- 
sible to  model  all  the  altimeter  errors  for  this  case  as  a 
single  bias  error. 

Attitude  Readout  Errors  — Since  scale  factor  errors 
are  one  of  the  dominant  sources  of  ARO  error  in  the  MESG 
(Ref.  6),  it  was  decided  to  characterize  the  ’ esidual  read- 
out errors  using  this  form.  As  a result  of  this  model  selec- 
tion, the  ARO  error  in  the  M-frame,  6y^  , that  appears  in 
Eq . (3.2-11)  was  expressed  as 


0 y±  ; i-1,2 


(3.3-2) 


Because  the  ARO  errors  propagate  into  MICRON  velocity  errors 
as  functions  of  the  RSA  direction  in  the  case  [see  Eq . (3.2-13)]  , 
the  states  in  the  error  analysis  program  that  model  these 
errors  are  products  of  the  terms  in  Eq.  (3.2-2)  and  the  com- 
ponents of  the  ^-vectors.  (Appendix  B enumerates  these  states 
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and  their  associated  dynamics.)  The  initialization  portion 
of  the  computer  program  accepts,  as  input  parameters,  the  co 
variance  statistics  of  the  6C^j's  in  Eq.  (3.3-2)  and  con- 
structs the  appropriate  states  required  to  implement  the 
effects  of  ARO  errors  as  described  in  Section  3.2. 


Drift  Rate  Errors  - The  model  for  residual  drift 
rate  of  the  MESG  spin  axis  includes  miscalibration  of  the 
case-fixed  bias  and  axial  mass  unbalance  drift  rate  coeffi- 
cients (see  Ref.  3 for  a discussion  of  these  drift  rate  terms) 
and  random  drift  rate  effects.  An  additional  drift  rate  er- 
ror that  results  from  computing  the  axial  mass  unbalance 
drift  rate  compensation  using  the  sensed  acceleration  is 
also  modeled.  (Since  the  sensed  acceleration  is  corrupted 
by  accelerometer  errors,  using  the  sensed  acceleration  in  the 
drift  rate  compensation  introduces  residual  drift  rate  errors 
which  drive  the  MICRON  navigation  errors.) 


The  parameters  which  describe  the  statistics  of  the 
case-fixed  bias  drift  rates  are  entered  into  the  error  analy- 
sis program  in  the  M-frame.  The  initialization  module  of 
the  program  then  computes  the  appropriate  covariance  matrix 
in  the  L-frame  to  model  these  errors.  To  simplify  the  mod- 
els, it  was  assumed  that  the  axial  mass  unbalance  drift  rate 
coefficient  and  its  calibration  error  were  the  same  for  both 
MESG's.  The  random  drift  rate,  which  was  modeled  as  a Markov 
process,  was  entered  into  the  program  in  the  L-frame.  It  is 
possible  to  summarize  the  model  for  residual  drift  rate  as 


L 

C 


♦ «Pg  £L  + Pg  1L 


(3.3-3) 
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where  e„,  and  eno  are  the  residual  case-fixed  drift  rate 

vectors  of  the  primary  and  secondary  MESG's.  Pg  is  the 

axial  mass  unbalance  coefficient  and  6P  is  the  calibration 

L ^ 

error  in  this  coefficient.  is  the  Markov  drift  rate 

vector . 


3.3.2  Program  Description 

The  structure  of  the  error  analysis  computer  program 
is  quite  simple  as  can  be  seen  from  the  flowchart  which  ap- 
pears in  Fig.  3.3-1.  The  modules  which  compose  the  program 
perform  initialization,  transition  matrix  computation,  up- 
dates to  model  a change  in  MICRON  attitude,  covariance 
propagation,  and  data  output.  The  initialization  module  ac- 
cepts the  input  data  which  describes  the  MICRON  error  par- 
ameters and  constructs  an  initia?  covariance  matrix,  PQ, 
to  represent  the  initial  conditions  on  the  system  state 
vector.  Also,  the  system  dynamics  matrix,  F,  and  the  con- 
tinuous noise  matrix,  Q,  are  assembled;  and,  data  defining 
the  initial  navigation  attitude,  program  time  step,  final 
time,  and  discrete  attitude  changes  are  input  and  stored  for 
use  by  other  modules. 

The  two  central  modules  of  the  error  analysis  pro- 
gram are  those  which  compute  the  transition  matrix  and  the 
discrete  noise  matrix  and  which  propagate  the  covariance 
matrix.  The  equations  which  accomplish  these  operations  are 
well-known  (Ref.  11),  but  will  be  repeated  here  for  complete- 
ness. The  transition  matrix,  $ , is  computed  as: 

$ = eFAt  (3.3-4) 

where  At  is  the  time  step  size  for  the  covariance  propaga- 
tion. The  discrete  noise  matrix,  Qk,  is  defined  by: 
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Qk  = J 4>  Q $T  dt  (3.3-5) 

0 

and  the  covariance  matrix  is  propagated  from  tfc  to  Vi 
(tk+l  - tk  + At)  by: 

P,  = ft  P.  <1>T  + Qk  (3.3-6) 

k+1  k K 

The  MICRON  error  analysis  computer  program  assumes 
that  the  MICRON  system  is  operating  in  a stationary  labora- 
tory environment,  except  for  discrete  changes  in  attitude. 

The  time  at  which  an  attitude  change  is  to  occur  and  the 
parameters  that  describe  the  attitude  change  are  program 
inputs.  The  MICRON  error  propagation  is  dependent  on  the 
orientation  of  the  RSA's  relative  to  the  gyro  case,  and  con- 
sequently, a change  in  MICRON  attitude  alters  the  driving 
error  sources.  The  attitude  change  module  properly  updates 
the  system  covariance  matrix  and  recomputes  the  dynamics 
and  noise  matrices  to  model  the  effects  of  an  attitude  change 
on  the  MICRON  navigation  errors. 

The  final  module  of  the  error  analysis  program  is  an 
output  module.  It  prints  the  diagonal  elements  of  the  co- 
variance  matrix  at  each  time  step  and  stores  those  data  points 
that  are  to  be  plotted.  Typically,  plots  have  been  generated 
for  the  time  history  of  the  rms  position,  velocity,  and  4>- 
angle  errors  of  the  MICRON  system  in  response  to  the  speci- 
fied instrument  error  sources.  The  next  section  will  present 
the  results  of  several  of  the  checkout  runs  that  were  made  to 
exercise  the  error  analysis  program. 

3.3.3  Covariance  Error  Analysis  Results 

The  error  analysis  computer  program  described  in  the 
previous  section  has  been  coded  and  exercised  to  verify  its 
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proper  operation.  Program  checkout  was  accomplished  by  set- 
ting all  error  sources  and  initial  condition  errors  to  zero 
except  one  and  ascertaining  that  the  single  error  source  in- 
duced the  appropriate  error  behavior.  Those  errors  which 
were  investigated  during  program  checkout  included  initial 
position  errors,  initial  velocity  errors,  initial  tji-angle 
errors,  bias  east  gyro  drift,  and  bias  north— axis  acceler- 
ometer error. 

At  the  conclusion  of  program  checkout,  a brief  in- 
vestigation was  conducted  using  the  error  analysis  program 
to  determine  the  navigation  errors  induced  by  those  error 
sources  which  are  unique  to  the  MICRON  system.  In  particular, 
the  effects  of  case-fixed  gyro  drift  and  scale  factor  ARO 
errors  were  studied  for  the  MICRON  system  in  the  stationary 
laboratory  environment. 

For  these  studies,  the  RSA  of  the  primary  MESG  was 
pointing  up  along  the  vertical  and  the  RSA  of  the  secondary 
MESG  was  pointing  north.  Thus,  the  RSA  vectors  in  L-frame 
coordinates  are: 


(3.3-7) 


The  attitude  of  the  MICRON  system  was  selected  so  that  the 
center  of  the  No.  1 electrode  was  parallel  to  the  earth's 
axis  of  rotation  (at  the  latitude  of  Anaheim,  California) 
and  that  the  x-axis  of  the  M-frame  was  in  the  plane  defined 
by  the  North  and  Down  axes  of  the  L-frame.  Consequently, 
the  transformation  to  the  L-frame  from  the  M-frame  is: 


3-26 


0 

0 

D 


U 

0 

u 

[1 

0 

c 

IJ 

0 

0 

0 

D 

0 

0 


0.816  -0.408  -0.408 

0 -0.707  0.707 

-0.577  -0.577  -0.577 


(3.3-8) 


Figure  3.3-2  presents  the  time  history  (for  a 12-hr 
period)  of  the  ^-angles  and  the  position  and  velocity  errors 
for  MICRON  when  the  only  error  is  a 0.005  deg/hr  drift  rate 
of  the  primary  MESG  about  the  x-case-fixed  axis.  Using  Eq. 
(3.3-8),  it  can  be  seen  that  a drift  rate  about  the  x-case- 
fixed  axis  should  cause  drifts  about  the  North  and  Down 
axes  in  the  L-frame.  The  results  in  Fig.  3.3-2  show  a 
ramping  growth  rate  for  *N  and  ^ which  is  the  behavior 
that  is  expected  due  to  constant  gyro  drift  rates.  The 
drift  rates  (through  the  interactions  of  the  ^-angle  mis- 
alignments with  the  resolution  of  gravity)  also  lead  to  pos- 
ition and  velocity  errors  which  exhibit  the  anticipated  be- 
havior including  Schuler  and  Foucault  oscillations  (Ref.  19) 
The  apparent  42  min  oscillation  of  the  North  velocity  error 
is  due  to  rectification  of  the  negative  portion  of  the 
Schuler  oscillation  occurring  because  covariance  results 
are  rms  (and  thus  strictly  positive)  numbers. 

The  navigation  errors  induced  by  an  x-axis  scale 
factor  ARO  error  for  the  primary  MESG  are  presented  in  Fig. 
3.3-3.  As  was  derived  in  Section  3.2  and  can  be  seen  from 
the  results,  ARO  errors  do  not  drive  the  ^-angle  equations, 
rather,  they  cause  velocity  errors  in  the  same  manner  as 
accelerometer  misalignments.  In  a similar  manner  to  the 
drift  rate  results,  the  apparent  42  min  oscillation  in  the 
position  and  velocity  errors  is  the  classical  Schuler  dy- 
namic behavior  masked  by  the  rectification  due  to  the  fact 
that  the  results  are  from  covariance  studies. 
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An  exhaustive  investigation  of  MICRON  system  perform- 
ance using  the  error  analysis  program  was  not  conducted  dur- 
ing the  present  effort  because,  as  was  discussed  previously, 
sufficient  data  has  not  been  collected  to  adequately  deter- 
mine statistical  models  that  describe  the  residual  instru- 
ment errors.  Rather,  a brief  investigation  to  determine 
the  error  propagation  characteristics  of  the  unique  MICRON 
error  sources  was  conducted.  The  insights  gained  from  this 
study  should  be  quite  useful  in  qualitative  evaluations  of 
the  test  data  collected  on  the  N57A  system  (a  developmental 
navigator,  using  two  MESG's  and  three  electro-magnetic  ac- 
celerometers, whose  purpose  is  to  demonstrate  the  strapdown 
MICRON  concept ) . 


The  MICRON  analysis  studies,  particularly  those  de- 
scribed in  this  chapter,  are  part  of  a continuing  TASC  ef- 
fort t"  support  the  MICRON  development  program.  It  has 
been  suggested  (at  the  beginning  of  this  chapter)  that, 
during  Phase  2 of  MICRON  development,  a sensor  data  bank 
be  established  to  develop  valid  statistical  models  of  re- 
sidual MICRON  instrument  errors  for  use  in  performance 
evaluation  and  aided-inertial  navigation  studies.  The  er- 
ror analysis  program,  described  in  this  chapter,  will  be  a 
very  important  tool  in  the  iterative  process  of  determining 
sensor  error  models  from  instrument  and  system  test  data. 


3-32 


TEST  SUPPORT 


In  addition  to  exercising  the  MCSP  and  developing  an 
error  analysis  program  to  statistically  investigate  MICRON 
performance,  analytic  support  was  provided  to  the  N57A  test 
program.  The  purpose  of  developing  the  N57A,  which  uses 
two  MESG's  as  its  orientation  reference  and  three  electro- 
magnetic accelerometers  (EMA's)  to  measure  specific  force, 
was  to  demonstrate  the  performance  attainable  with  a strap- 
down  MESG  inertial  navigator.  The  N57A  was  not  designed  to 
demonstrate  MICRON'S  ultimate  low  cost  and  small  size. 


Section  4.1  will  review  the  N57A  Test  Plan  and  dis- 
cuss recommended  improvements  which  resulted  from  this  review. 
The  results  of  the  N57A  tests  conducted  by  Autonetics  will 
be  presented  in  Section  4.2  and  observed  performance  will 
be  compared  with  that  delineated  in  the  N57A  specification. 
Section  4.3  will  summarize  the  results  and  identify  areas 
requiring  additional  development. 


4.1  REVIEW  OF  THE  N57A  TEST  PLAN 


The  goal  of  the  N57A  test  program  was  to  determine 
the  system  performance  under  both  ambient  and  environmental 
conditions.  Table  4.1-1  presents  the  performance  targets  for 
the  N57A  while  Table  4.1-2  describes  the  environmental  con- 
ditions which  the  system  must  withstand  (Ref.  6).  In  addi- 
tion to  system  level  tests  of  the  N57A,  the  Test  Plan  (Ref.  20) 
describes  a complete  sequence  of  module,  subassembly  and  sub- 
system tests.  The  purpose  of  this  testing  was  to  verify  the 
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TABLE  4,1-1 

PERFORMANCE  TARGETS  FOR  THE  N57A 


Horizontal  Position  Error 

Horizontal  Velocity  Error 
Vertical  Velocity  Error 
Azimuth  Error 
Roll,  Pitch  Errors 
Tilt  Error 


1 nm/hr,  CEP  (Least  Squares, 

straight  line  fit) 

5 ft/sec,  time  rms,  per  axis 

2 ft/sec,  lo 

4 min,  1 a 

4 min,  la 

1 min,  time  rms,  per  axis 
1 hr  from  0°F 
16  hrs 
30  days 


Reaction  Time 
Calibration  Time 
Calibration  Stability 


integrity  of  the  component  functions  at  an  intermediate  level 
to  reduce  the  problems  associated  with  system  integration. 

The  N57A  system  level  tests  begin  during  the  last 
phase  of  system  integration  and  continue  through  a lab  test 
and  a van  test.  Flight  tests  of  the  N57A,  conducted  by  both 
Autonetics  and  the  Central  Inertial  Guidance  Test  Facility 
(CIGTF)  at  Holloman  AFB,  are  planned,  following  the  van  test 
program.  Briefly  summarized,  the  tests  to  be  conducted  dur- 
ing each  phase  are: 

• System  Integration  Tests 

Align  the  N57A  at  a minimum  of  three  headings 
and  attitudes. 

Perform  4-hr  navigation  runs  at  a minimum  of 
three  different  headings,  changing  heading 
and  IMU  tilt  after  two  hrs  during  some  of  the 
runs. 

• Lab  Tests 

Four-hr  navigation  under  Scorsby  motion. 

The  N57A  will  be  exposed  to  environments 
including  linear  vibration,  correlated  angu- 
lar and  linear  vibration,  shock,  and  low 
temperature. 

• Van  Tests 

Align  in  the  van  at  several  headings. 

Perform  4-hr  navigation  runs  using  both 
triangular  and  constant  heading  courses. 

The  system  integration  and  lab  test  phases  of  the  test 
program  are  comprehensive  and  well-defined.  In  particular,  the 
procedure  which  requires  shakedown  tests  in  each  different  en- 
vironment, where  the  gyros  are  suspended  but  not  spinning,  prior 
to  navigation  tests  in  the  environment  is  most  prudent.  This 
procedure  will  hopefully  minimize  the  danger  of  secondary  fail- 
ure of  the  MESG.  Successful  completion  of  this  test  phase 
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should  accomplish  the  dual  objectives  of  evaluating  N57A  per- 
formance under  the  required  conditions  and  qualifying  the  sys- 
tem for  van  testing. 

Upon  completing  the  review  of  the  proposed  lab  test 
program,  TASC  made  several  suggestions  concerning  potential 
improvements  to  the  proposed  test  sequence.  Although  Autonetics 
proposed  a minimum  testing  philosophy  with  no  redundant  or  un- 
necessary tests,  it  was  suggested  that  static  navigation  runs 
be  scheduled  between  each  different  set  of  environmental  tests. 
The  reason  for  these  additional  baseline  navigation  tests  was 
to  verify  that  the  system  was  moved  to  the  new  test  equipment 
properly  and  to  ascertain  that  no  permanent  performance  degrada- 
tion had  been  induced  by  the  previous  environmental  test.  These 

tests  will  also  generate  more  performance  data  on  the  N57A  sys- 
tem. 


The  longest  navigation  test  run  that  had  been  proposed 
by  Autonetics  was  four  hrs  in  duration.  TASC  suggested  that  at 
least  one  longer  navigation  test  lasting  6-8  hrs  be  conducted 
during  the  lab  tests.  This  longer  duration  navigation  run 
should  be  used  to  demonstrate  that  there  are  no  error  sources 
peculiar  to  a MESG  navigator  that  appear  only  after  extended 
operation.  (Eight  hours  was  selected  for  the  maximum  duration 
of  the  navigation  testing  since  that  time  interval  will  encom- 
pass most  all  mission  applications  for  a moderate  accuracy  air- 
craft navigator).  Also,  a navigation  test  of  this  duration  will 
allow  preliminary  investigation  of  the  periodic  long  term  error 
behavior  associated  with  unaided  inertial  navigation  (see  Ref. 

15  for  a discussion  of  inertial  system  error  dynamics). 

The  Test  Plan  published  by  Autonetics  did  not  describe 
in  great  detail  the  van  tests  that  were  to  be  performed. 
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Following  review  of  the  proposed  van  tests,  TASC  suggested 
several  improvements  to  the  planned  test  procedures  These 
suggestions  include: 

• Conducting  long  duration  (6  to  8 hr)  van 
runs  . 

• Implementing  an  altimeter  and  an  active 
vertical  channel. 

• Providing  a velocity  reference  against  which 
to  evaluate  N57A  performance  while  the  van 
is  in  motion. 

• Defining  a test  course  which  has  a significant 
variation  in  altitude  (i.e.  , driving  the  test 
van  into  mountainous  areas). 

The  reason  for  suggesting  longer  duration  van  runs 
are  the  same  as  those  discussed  for  longer  navigation  tests  in 
the  lab.  Including  an  altimeter  and  an  active  vertical  channel 
in  the  van  will  allow  checkout  of  these  systems  prior  to  moving 
the  N57A  into  the  flight  test  aircrafts.  Providing  a velocity 
reference  permits  velocity  performance  evaluation  while  the  van 
is  in  motion  which  may  facilitate  identification  of  the  causes 
of  N57A  navigation  errors.  Using  a test  course  with  significant 
altitude  variations  will  enhance  the  van  test  program  because  it 
allows  adequate  exercising  of  the  vertical  channel  and  it  induces 
an  attitude  change  of  the  N57A  during  the  steep  climbs. 

In  summary,  TASC  has  reviewed  the  N57A  Test  Plan  and 
made  comments  and  suggestions  in  areas  of  potential  improve- 
ment . Many  of  the  suggestions  have  been  incorporated  into 
the  testing  of  the  N57A  conducted  by  Autonetics.  In  particu- 
lar, all  the  suggested  modifications  of  the  van  test  program 
were  to  be  incorporated  into  the  second  sequence  of  van  tests 
conducted  in  March  1974. 
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4.2  N57A  TEST  RESULTS 

The  lab  and  van  tests  of  the  N57A  system  included  90 
navigation  runs  under  a large  number  of  different  conditions. 
The  results  of  these  tests  presenting  time  history  plots  of 
position  and  velocity  errors  are  contained  in  Refs.  21  and 
22.  Rather  than  reproducing  all  these  results  here,  Table 
4.2-1  presents  a summary  of  the  number  of  tests  of  each  type 
that  were  conducted.  This  section  will  briefly  discuss  the 
results  of  each  different  environmental  test  and  then  com- 
ment on  the  performance  of  the  N57A  system  during  the  entire 
test  program. 


TABLE  4.2-1 

SUMMARY  OF  N57A  LAB  AND  VAN  TESTS 


Type  of  Test 

Number  of 
Tests 

Laboratory  — Stationary 

20 

Laboratory  — Stationary  with  Discrete 
Attitude  Changes 

12 

Scorsby  Motion 

4 

Linear  Vibration 

9* 

Correlated  Linear  and  Angular  Vibration 

3 

Angular  Rates 

5 

Shock 

4** 

Low  Temperature 

1 

Van  — Stationary 

19 

Van  — Moving 

17 

♦ 

Four  tests  were  conducted  with  the  shaker  field  on, 
+)(tbut  with  no  vibration  (see  text) 

System  was  not  operating  during  shock 


0 


j, 

[ 

[ 

I 

f 

l 

1 

0 

0 

l 
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4.2.1  Environmental  Test  Results 


Since  the  performance  goal^s  are  specified  as  en- 
semble statistics,  relating  N57A  performance  during  the  en- 
vironmental tests  to  the  performance  targets  presented  in 
Table  4.1-1  is  difficult  because  only  a few  tests  were  con- 
ducted in  each  environment.  In  a single  test  it  is  always 
possible  to  get  a performance  sample  which  is  significantly 
worse  than  true  average  performance  of  the  system.  Only 
when  a large  number  of  samples  (tests)  are  taken,  can  the 
results  be  viewed  with  confidence  as  a true  indication  of 
system  performance.  Because  of  this  small  sample  size  prob- 
lem, the  N57A  environmental  test  results  are  discussed  only 
qualitatively,  pointing  out  potential  performance  problems 
that  appear  to  be  present  in  the  few  tests  performed. 

The  results  of  the  four  Scorsby  tests  of  the  N57A 
that  consist  of  four  hours  of  navigation  mode  operation  in  the 
presence  of  continuous  Scorsby  motion  indicate  that  this  en- 
vironment does  not  induce  any  significant  errors.  The  per- 
formance of  the  N57A  during  the  Scorsby  testing  is  quite 
similar  to  performance  obtained  during  stationary  lab  and 
van  tests.  Also,  all  four  navigation  runs  during  Scorsby  mo- 
tion have  similar  error  magnitudes  indicating  that  these 
samples  are  likely  representative  of  the  N57A  performance 
during  this  motion  and  not  anomalous  results  which  can  always 
occur  in  single  trials  with  systems  driven  by  stochastic  error 
sources. 


The  performance  of  the  N57A  during  the  linear  vibra- 
tion tests  was  not  as  ^ood  as  during  the  other  environments. 
These  vibrations,  which  consisted  of  sweeps  between  2 and  10 
g's  in  amplitude  as  the  frequency  was  varied  from  5 to  2000  Hz, 
caused  particularly  poor  performance  when  they  were  applied 
along  the  vertical  axis.  Investigation  of  the  test  set  up 
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led  to  the  discovery  that  the  shaker,  which  generated  the  vi- 
bration environment,  produced  a sizable  magnetic  field.  This 
magnetic  field  caused  a noticeable  speed  reduction  of  the  MESG 
rotors  when  the  vibration  was  applied  along  the  vertical  axis. 

Attempting  to  separate  the  effects  of  the  shaker's 
magnetic  field  from  the  vibration  environment,  navigation  runs 
with  the  shaker  field  on  but  no  vibrations  were  made  in  each 
test  configuration.  (Tests  were  made  which  applied  vibra- 
tions along  each  of  the  vertical,  lateral  horizontal,  and 
longitudinal  horizontal  axes.)  These  tests  indicated  that 
performance  degradation  occurred  with  and  without  the  vibra- 
tion as  long  as  the  shaker  magnetic  field  was  present.  The 
performance  degradation  was  most  severe  in  the  vertical  vibra- 
tion test  set  up.  The  N57A  performance  during  the  horizontal 
vibration  test  configuration  did  not  appear  to  be  affected  very 
strongly  by  the  shaker's  field.  Summarizing  the  linear  vibra- 
tion tests,  N57A  performance  in  this  environment  has  not  been 
clearly  demonstrated  due  to  disturbances  induced  by  the  test 
procedure,  but  it  seems  likely  that  a disturbance-free  test 
environment  would  demonstrate  satisfactory  navigation  perform- 
ance in  the  presence  of  linear  vibrations. 

Neither  the  correlated  linear  and  angular  vibration 
tests,  nor  the  shock  tests  identified  any  significant  weakness 
in  the  N57A  system  performance.  The  results  of  correlated 
motion  tests,  which  consisted  of  one  test  at  each  of  three  fre- 
quencies (64  Hz,  6.4  Hz,  and  the  resonant  frequency  of  the 
shock  mounts),  were  all  consistent  with  N57A  performance  goals. 
Although  not  unacceptably  large,  the  worst  errors  (radial  error 
rate  = 1.49  nm/hr,  rms  velocity  error  ■ 2.995  fps)  occurred 
during  64  Hz  correlated  motion.  This  frequency  is  the  solution 
rate  of  the  N57A  fast  cycle  equations  (Ref.  14)  during  which 
the  accelerometers  are  sampled  and  projected  on  the  inertial 
frame  defined  by  the  MESG  spin  axis  attitude  readouts.  It  would 


be  expected  that  correlated  input  motion  at  this  frequency  would 

be  a source  o,  error.  The  shoe*  input  of  a 1*.  » — « 1 

sine  wave  was  applied  while  the  N57A  was  operating  but  without 
the  MESG  rotors  spinning.  This  environment  did  not  indu e any 
system  failures  and  no  performance  degradation  was  detected 
lowing  the  exposure  to  the  shock  environment . 


The  performance  of  the  N57A  at  low  temperature  (0  F)  was 
comparable  to  the  other  stationary  lab  test  resuits.  Theje^ 
for  the  low  temperature  test  was  to  determine  t e Ie 
o,  the  system  at  oV  For  this  test  the  reaction  time  w s 1 4 
min  (Ref.  23)  which  is  twice  the  target  of  ~ 

Table  4.2-1.  Reaction  time  is  an  area  requ  r ng 
improvement  from  the  N57A  to  the  MICRON  system  and  is  a tas 
receiving  considerable  research  and  developmen e”"V  and 
Autonetics  (Ref.  24).  Incorporation  of  automatic  « 
heating  should  Improve  the  system  reaction  tlme  b* 
the  slow  processes  of  manual  polhode  damping  and  rotor 


The  angular  rate  testing  included  one  run  with  three 
hours  of  continuous  12°/sec  yaw  rate  and  three  one-hour  tests 
with  short  duration  (only  long  enough  to  get  the  test  table  up 
to  the  desired  speed  and  stopped  aRain)  high  rate  inputs  a 


to  the  aesirea  ~ °,  in 

30  min.  The  high  rate  inputs  were  40«  sec  n yaw  0 e n 


roll  and  100°/sec  in  pitch.  The  angular  rate  inputs  did  not 
appear  to  cause  any  significant  degradation  o,  N57A  navigation 
accuracy.  In  summary,  the  N57A  accuracy  does  not  appea 
reduced  by  Scorsby  motion,  correlated  linear  and  angular  vibr  - 
tion,  shock,  low  temperature  or  high  angular  ra  • 
served  performance  degradation  in  the  presence  of  linear  v 
tion  is  likely  due  to  the  test  set  up  and  not  to  the  vibration 
“ tut  a modified  test  configuration  would  be  needed  to  verify 
the  vibration  performance.  The  reaction  time  of  the  N57A  ron. 
low  temperature  did  not  meet  the  goal,  indicating  t a a 
tional  effort  is  required  in  this  area. 


I 
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4.2.2  N57A  Performance  Evaluation 


To  quantitatively  estimate  the  N57A  system  performance 
from  the  lab  and  van  test  results,  all  90  navigation  runs  were 
considered  to  provide  a data  base  large  enough  to  be  statisti- 
cally significant.  Figures  4.2-1  and  4.2-2  summarize  on  his- 
tograms the  radial  error  rate  and  time  rms  velocity  error  re- 
sults observed  by  Autonetics  during  the  test  program.  Different 
symbols  are  used  in  the  figures  to  distinguish  stationary  lab 
tests,  lab  environmental  tests,  stationary  van  tests  and  mov- 
ing van  tests. 

The  horizontal  position  error  growth  rate  performance 
target  for  the  N57A  is  1 nm/hr  CEP  (least  squares  straight 
line  fit  to  radial  error)  which  means  that,  on  the  average, 

50%  of  the  time  the  observed  radial  error  rate  will  be  less 
than  1 nm/hr.  The  estimate  of  the  position  error  performance 
of  the  N57A  was  computed  by  tabulating  the  computed  slopes  of 
a least  squares  straight  line  fit  to  the  radial  position  error 
time  histories  for  the  test  results  and  selecting  the  median 
value.  This  method  yielded  a value  of  0.87  nm/hr  as  the  CEP 
rate  for  the  N57A. 

The  N57A  performance  target  for  velocity  accuracy  is 
a lo  value  of  5 ft/sec,  time  rms  per  axis.  A lo  performance 
number  implies  that,  on  the  average,  68%  of  the  time  the  sys- 
tem accuracy  will  equal  or  better  the  specified  value.  The 
68  percentile  velocity  error  performance  was  computed,  in  a 
similar  manner  to  the  position  computation,  to  be  3.15  ft/sec. 

Measurement  of  the  attitude  errors  of  the  N57A  during 
a navigation  run,  particularly  when  motion  is  involved,  is 
difficult  to  accomplish  due  to  problems  in  defining  a refer- 
ence against  which  N57A  outputs  can  be  compared.  Consequently, 
very  little  attitude  accuracy  data  was  collected  during  the 
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RADIAL  ERROR  RATE  («m/hr) 


Figure  4.2-1 


N57A  Position  Error  Histogram 
(Data  Provided  by  Autonetics) 
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Figure  4.2-2 


N57A  Velocity  Error  Histogram 
(Data  Provided  by  Autonetics) 
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test  program.  For  about  one-third  of  the  tests,  the  initial 
azimuth  error,  which  is  indicative  of  system  performance 
during  alignment,  was  measured.  Figure  4.2-3  presents  these 
results  in  histogram  form.  (Different  symbols  are  used  to 
distinguish  lab  test  from  van  test  results.  No  distinction 
is  made  between  moving  and  stationary  test  results  since  the 
N57A  is  stationary  during  the  alignment  phase  of  all  tests.) 
Computing  the  la  azimuth  error  from  the  test  data  as  the  68 
percentile  yields  an  estimate  of  initial  azimuth  accuracy  for 
the  N57A  of  6.1  min.  The  performance  target  for  azimuth  ac- 
curacy was  4 min,  la  (see  Table  4.1-1). 
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Figure  4.2-3 


N57A  Initial  Azimuth  Error  Histogram 


No  test  results  were  collected  on  the  accuracy  of  the 
vertical  channel  during  the  N57A  lab  and  van  tests  because 
the  system  was  not  integrated  with  an  altimeter  during  this 
phase  of  the  test  program.  Rather,  the  input  to  the  vertical 
channel  was  pinned  at  a constant  value.  During  the  second 
set  of  N57A  van  tests,  an  altimeter  was  scheduled  to  be  in- 
cluded and  the  vertical  performance  of  the  N57A  was  to  be 
evaluated.  Prior  to  the  beginning  of  these  tests,  the  N57A 
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was  modified  to  incorporate  completely  automatic  startup  cap- 
ability. This  modification  induced  a hardware  problem  (the 
MESG  feed-through  pins  became  magnetized  during  startup) 
which  prevented  any  useful  performance  data  from  being  col- 
lected during  these  van  tests.  It  was  possible  to  check  out 
the  altimeter  interface,  but  evaluation  of  the  vertical  chan- 
nel performance  was  postponed  until  the  flight  tests. 

Evaluation  of  the  stability  of  the  N57A  calibration 
parameters  during  the  test  program  is  difficult  because  hard- 
ware problems  prevented  extended  operation  of  the  system  be- 
tween calibrations.  Table  4.2-2  presents  the  calibration 
history  of  the  N57A  including  the  reason  for  each  recalibra- 
tion and  the  number  of  tests  performed  in  the  interval.  The 
heading  sensitivity  which  precipitated  the  recalibration  fol- 
lowing Test  No.  29  was  isolated  to  a calibration  parameter 


shift.  All  the  other  recalibrations  made  during  the  lab  and 
van  testing  were  due  to  causes  other  than  calibration  param- 
eter shifts.  Thus,  the  one  sample  that  is  indicative  of 
N57A  calibration  stability  is  28  days.  (The  target  for  cal- 
ibration stability  was  30  days.)  Additional  testing  will 
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TABLE  4.2-2 

CALIBRATION  STABILITY  RESULTS 


CALIBRATION  DATS 
(Approx la* t.) 

SYSTEM  OPERATION 

Dates 

Test  Numbers 

8/10/73 

9/10/73  - 9/  7/73 

1-29 

9/14/73 

9/14/73  - 9/25/73 

30-47K 

10/24/73 

10/24/73  - 10/26/73 

48-52 

11/  9/73 

11/  9/73  - 11/27/73 

53-71 

1/24/74 

1/24/74  - 1/23/74 

72-80 

REASON  FOR 
RECALIBRATION 


TIME  BETWEEN 
CALIBRATIONS 
(Days) 


System  developed 
heading  censltlvlty 

Failed  hardware  In 

system 

System  failed  follow- 
ing cold  soak 

Test  sequence  completed 


Automatic  startup  prob- 
lems 


28 

12 

2 

18 

1 


ti 
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clearly  be  required  to  develop  confidence  that  the  stability 
of  the  MESG  compensation  parameters  is  adequate  to  achieve 
the  desired  calibration  intervals. 
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4.3  TEST  SUPPORT  SUMMARY  AND  CONCLUSIONS 

In  support  of  the  N57A  test  program,  TASC  has  reviewed 
the  proposed  system  level  lab  and  van  test  plans  and  made  rec- 
ommendations concerning  potential  improvements  to  the  planned 
tests  and  procedures.  Several  of  these  suggestions  which  in- 
cluded baseline  navigation  runs  between  environmental  tests, 
longer  duration  navigation  runs,  and  incorporation  of  an  al- 
timeter and  a velocity  reference  in  the  mobile  test  vehicle 
were  adopted  by  Autonetics  for  portions  of  the  test  program. 

The  results  of  the  N57A  lab  and  van  tests  were  evalu- 
ated to  determine  the  system  performance.  This  performance 
may  be  summarized  as 

• Radial  Error  Rate  CEP  = 0.87  nm/hr 

• Horizontal  Velocity  Error,  lo  =3.2  ft/sec 

• Initial  Azimuth  Error,  lo  = 6.1  min 

• Reaction  Time  from  0°F  = 124  min 

• Calibration  Stability  £ 28  days 
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The  N57A  system  accuracy  during  the  test  program  compares 
reasonably  well  with  the  performance  targets  presented  in 
Table  4.1-1.  Those  areas  which  the  test  program  has  identi- 
fied as  needing  additional  testing  and  development  are  reac- 
tion time  and  calibration  stability.  Also,  flight  test  per- 
formance still  must  be  demonstrated.  There  are  some  system 
level  error  sources  which  will  be  more  visibly  excited  by 
the  high  acceleration  and  velocity  environments  of  an  air- 
craft [in  particular,  the  (£  x £)  term  in  Eq.  (3.2-2)]. 
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Thus,  it  may  be  anticipated  that  larger  navigation  errors 
will  be  observed  during  the  flight  tests  than  were  seen  dur- 
ing lab  and  van  testing. 

Comparing  the  observed  N57A  performance  with  the  error 
budget  developed  for  this  system  early  in  the  MICRON  Program 
points  out  a conflicting  set  of  circumstances.  While  the 
N57A  navigation  accuracy  was  generally  better  than  the  tar- 
get, the  MESG  drift  rate  and  attitude  readout  residuals  were 
larger  than  allowed  by  the  error  budget.  This  discrepancy 
is  likely  due  to  the  use  of  improper  statistical  models  for 
the  MESG  error  residuals  in  generating  the  error  budget.  Be- 
cause of  the  high  order  of  the  MESG  compensation  models,  the 
error  residuals  are  likely  of  high  harmonic  order  in  case- 
rotor  orientation.  Consequently,  the  vibration  environment 
to  which  a strapdown  navigator  is  exposed  may  reduce  the 
effective  residual  errors  by  "dithering"  the  RSA  vector  in 
case  coordinates.  It  is  suggested  that  future  studies  be 
undertaken  to  attempt  to  determine  valid  statistical  models 
for  the  MESG  error  residuals.  These  models  are  required 
not  only  for  reconciliation  of  the  MICRON  error  budget  but 
also  for  development  of  Kalman  filtering  techniques  for  in- 
tegrating MICRON  with  various  navaids  such  as  doppler  radar, 
OMEGA,  and  navigation  satellites. 

Summarizing,  the  N57A  lab  and  van  tests  have  accom- 
plished the  goal  of  demonstrating  1 nm/hr  navigation  accuracy 
using  MESG's  in  a strapdown  mode.  The  environmental  testing 
conducted  tends  to  indicate  the  navigation  accuracy  is  not 
significantly  degraded  by  exposure  to  the  test  environments. 
(Linear  vibration  resistance  requires  further  demonstration 
due  to  test  eauipment  induced  magnetic  fields.)  Further 
testing  is  planned  to  demonstrate  the  performance  of  the 
N57A  in  both  fighter  and  transport  aircraft  and  helicopters. 
Based  on  the  results  of  the  lab  and  van  tests,  the  two-hour 


reaction  time  in  a low  temperature  environment  and  the  cali- 
bration stability  are  the  performance  areas  requiring  great- 
est improvement  to  develop  an  inertial  navigation  system 
suitable  for  operational  deployment. 
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5. 


SUMMARY  AND  CONCLUSIONS 


The  results  of  three  closely  related  system-level 
analyses  conducted  by  TASC  in  support  of  the  MICRON  develop- 
ment program  have  been  presented.  These  studies  consisted 

of : 

• Exercising  the  previously  developed  MICRON 
Computer  Simulation  Program  (MCSP)  to  eval- 
uate the  effects  of  MESG  instrument  varia- 
tions and  severe  environments  on  navigation 
performance . 

• Developing  models  to  describe  the  propaga- 
tion of  instrument  errors  into  navigation 
system  errors  and  generating  a covariance 
analysis  computer  program  to  conduct  sta- 
tistical performance  predictions  for  MICRON. 

• Supporting  the  test  program  of  the  N57A 
demonstration  system  by  reviewing  the  Test 
Plan  developed  by  Autonetics  and  by  evalu- 
ating the  laboratory  and  van  test  data. 

The  MCSP  studies  investigated  the  effects  on  MICRON 
performance  of  reducing  the  preload  charge  level  on  the 
MESG's,  eliminating  the  instrument  turn-on  repeatability 
errors,  and  subjecting  the  system  to  high  vibration  and 
roll  rate  environments.  A 15%  performance  improvement  was 
observed  in  the  simulation  results  when  the  preload  charge 
level  was  reduced  by  30%.  Comparable  performance  improve- 
ment was  also  predicted  by  the  MCSP  studies  when  only  the 
MESG  calibration  errors  (gyro  turn-on  repeatability  errors 
were  eliminated)  were  used  to  drive  the  system.  The  MCSP 
results  indicated  that  neither  high  vibration  nor  high  roll 
rate  and  acceleration  environments  should  prevent  the  MICRON 
system  from  achieving  its  specified  navigation  accuracy. 


5-1 


A covariance  analysis  computer  program  (to  be  used 
for  statistical  error  propagation  studies)  was  developed  to 
supplement  the  MCSP  and  to  increase  the  capability  to  simu- 
late MICRON  system  performance.  The  derivation  of  the  MICRON 
error  equations  has  been  presented,  the  structure  of  the  co- 
variance  analysis  program  has  been  described,  and  the  results 
of  a brief  investigation  to  determine  the  propagation  char- 
acteristics of  those  error  sources  unique  to  MICRON  have  been 
discussed.  Because  there  is  insufficient  data  available  to 
allow  development  of  valid  statistical  models  describing  re- 
sidual MESG  errors,  it  was  assumed  during  the  present  studies 
that  the  errors  were  due  to  miscalibration  of  several  of  the 
drift  rate  and  ARO  error  compensation  coefficients. 

It  is  strongly  suggested  that  development  of  models 
to  describe  the  residual  MESG  errors  be  given  high  priority 
during  future  MICRON  Analysis  studies  since  these  models  of 
the  MICRON  system  are  needed  for  use  in  modern  data  proces- 
sing algorithms  (e.g.,  Kalman  filtering)  which  will  combine 
MICRON  data  with  the  navaids  that  are  available  in  a modern 
avionics  suite.  To  aid  in  this  effort,  it  was  suggested 
that  a sensor  data  bank  in  which  to  store  data  collected  on 
all  MESG's  that  are  manufactured  be  established.  Residual 
error  model  determination  is  an  iterative  process  during 
which  test  data  is  analyzed,  models  are  postulated,  and 
model  predicted  system  performance  is  compared  with  test 
data.  The  covariance  analysis  program  developed  in  this  re- 
port will  be  a very  useful  tool  in  this  effort. 

In  support  of  the  N57A  test  program,  TASC  reviewed 
the  Test  Plan  document  and  evaluated  the  test  data  collected. 
The  review  of  the  Test  Plan  suggested  several  areas  where 
improvement  could  be  made.  In  particular,  with  respect  to 
the  van  test  program,  the  following  suggestions  were  made: 
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• Conduct  long  duration  (6  to  8 hr)  van 
runs . 

• Implement  an  altimeter  and  an  active 
vertical  channel. 

• Provide  a velocity  reference  in  the 
van  against  which  to  evaluate  N57A  per- 
formance . 

• Include  significant  altitude  variations 
in  the  van  test  courses. 

Many  of  these  suggestions  were  adopted  by  Autonetics  for  the 
second  set  of  van  tests  planned  for  the  N57A. 

Evaluation  of  the  data  collected  during  the  N57A 
laboratory,  environmental,  and  van  test  programs  indicates 
that  the  observed  system  performance  achieves  most  of  the 
targets  that  were  established  for  the  system.  The  computed 
system  accuracy  based  on  all  the  tests  performed  was 

• Radial  error  rate  CEP  = 0.87  nm/hr 

• Horizontal  velocity  error  (la)  = 3.2  ft/sec 

None  of  the  environments  to  which  the  N57A  was  subjected  iden- 
tified any  significant  performance  deficiency,  although  vi- 
bration performance  has  not  been  completely  verified  due  to 
test  equipment  problems.  Those  areas  which  have  been  iden- 
tified by  ttoe  results  of  the  test  program  as  requiring 
additional  testing  and  development  are  reaction  time  and 
calibration  stability. 

A reevaluation  of  the  MICRON  design  error  budget 
should  be  conducted  in  light  of  the  results  of  the  N57A  test 
program.  Due  to  the  lack  of  valid  models  to  describe  re- 
sidual MESG  errors,  the  design  error  budget  apparently  im- 
posed unnecessarily  severe  specifications  on  MESG  residual 
errors.  (Although  the  observed  system  performance  was  con- 
sistent with  the  targets,  the  residual  instrument  errors 
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were  generally  larger  than  the  values  allowed  in  the  error 
budget.)  Additional  analytic  effort  will  be  required  to 
resolve  this  issue  and  define  a revised  error  budget  for 
the  final  MICRON  design. 
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APPENDIX  A 


ERROR  EQUATIONS  FOR  MICRON'S 
VERTICAL  CHANNEL  DAMPING 


The  vertical  channel  damping  of  the  MICRON  system  is 
implemented  by  using  an  altimeter  to  generate  a vertical  posi 
tion  error  signal  which  is  fed  back  to  the  input  of  the  verti 
cal  velocity  and  position  integrators.  The  altimeter  output 
is  not  used  in  the  gravity  calculation.  Viewing  the  vertical 
channel  as  uncoupled,  Fig.  A-l  represents  the  mechanization 
of  the  vertical  channel  for  MICRON  (this  corresponds  to  the 
< as  o case  discussed  in  Ref.  26).  The  gains  and  kg  in 
Fig.  A-l  are,  for  the  MICRON  system, 

kx  = 2/t 

. o 2 . 1 

k2  = 2WS  t2  ^2 

where  t * 100  sec  and C * 0.7  (Ref.  21). 

As  derived  in  Ref.  26,  the  error  diagram  that  cor- 
responds to  the  MICRON  vertical  channel  mechanization  is 
presented  in  Fig.  A-2.  From  this  error  diagram,  it  is  pos- 
sible to  write  the  uncoupled  vertical  channel  error  equa- 
tions as 

6’rD  = 6vD  " kl(6rD  + <A_3) 

fiVn  = 6fn  + 2w_2  6rn  - k„( 6rn  + «hR) 


(A-l) 

(A-2) 


( A-4  ) 


GRAVITY 

CALCULATION 


f „ - VERTICAL  SPECIFIC  FORCE 

vr>  - VERTICAL  VELOCITY 

r„-  VERTICAL  POSITION  (POSITIVE  DOWN) 

hR-  ALTIMETER  OUTPUT  (POSITIVE  UP) 


Vertical  Channel  Mechanization 
for  the  MICRON  System 


Figure  A-l 


R-1«U2 


MICRON  Altitude  Channel  Error  Diagram 


Figure  A-2 


T 


The  altimeter  output  and  vertical  position  are  defined  with 
their  positive  directions  being  opposite  (altimeter  output 
is  positive  up,  while  vertical  position  is  positive  down  in 


a NED  frame)  which  accounts  for  the  sum  (6rD  + 6hR)  being  fed 


back  in  the  error  equations. 


The  term  6fD  in  Eq.  (A-4)  is  the  third  element  of 
6fL  defined  in  Eq.  (3.2-12)  and  includes  effects  of  both  the 
accelerometer  and  ARO  errors  that  excite  the  vertical  chan- 
nel. Equations  (A-3)  and  (A-4),  with  the  inclusion  of  the 
coupling  terms,  are  the  equations  that  model  the  vertical 
channel  in  the  error  analysis  computer  program.  Appendix  B 
contains  a complete  listing  of  the  state  vector  and  system 
dynamics  matrix  that  describe  the  MICRON  navigation  perform- 
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APPENDIX  B 

MODEL  DEFINITION  FOR 
ERROR  ANALYSIS  COMPUTER  PROGRAM 

This  appendix  presents  the  state  vector  and  dynamics 
matrix  implemented  in  the  covariance  error  analysis  computer 
program.  In  Chapter  3,  the  MICRON  error  equations  were  de- 
rived in  vector-matrix  form  and  the  models  describing  the 
residual  instrument  errors  that  drive  the  error  equations 
were  selected.  Table  3.3-1  indicates  that,  in  addition  to 
the  9 system  states,  33  error  sources  were  included  in  the 
MICRON  system  error  model.  A vector  of  45  states  is  required 
to  represent  the  MICRON  navigation  performance  in  a station- 
ary, laboratory  environment. 

Table  B-l  defines  the  elements  of  this  state  vector 
which  is  the  basis  of  the  error  analysis  computer  program. 
Since  it  was  decided  to  perform  the  error  analysis  in  the 
local-level,  north-east-down  L-frame,  all  the  elements  of 
the  system  state  vector  are  coordinatized  in  this  frame. 

The  describing  parameters  of  those  instrument  errors  which 
occur  in  the  gyro  case  fixed  M-frame  (case  fixed  bias  drift 
rates,  accelerometer  errors,  and  scale  factor  ARO  errors) 
are  input  relative  to  this  frame  and  transformed  by  the 
error  analysis  program  initialization  module  to  generate 
equivalent  errors  in  the  L-frame. 

For  the  stationary,  laboratory  environment,  the 
gravity  vector  is  always  along  the  down-axis  of  the  L-frame. 
This  fact  has  been  used  to  reduce  the  number  of  states  re- 
quired to  model  the  MICRON  system  (for  example,  only  the 
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TABLE  B-l 

STATE  VECTOR  FOR  ERROR  ANALYSIS  PROGRAM 


misalignments 


navigator  velocity  trror 


navigator  poaltlon  error 
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4C2D  y20 

*ie 

m 
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"vertical"  accelerometer  scale  factor  is  required).  Since 
drift  rate  and  ARO  errors  propagate  into  MICRON  navigation 
errors  as  functions  of  the  RSA  direction  vectors  (see  Sec- 
tion 3.2),  a large  number  of  states  are  required  to  properly 
model  these  effects.  Although  the  error  terms  in  these 
cases  (case  fixed  drift  rate  or  ARO  scale  factor  error) 
have  no  dynamics,  the  states  required  in  the  model  are  gov- 
erned by  the  dynamics  of  the  RSA  vectors.  For  the  present 
study,  the  motion  of  the  RSA  vectors  relative  to  the  L-frame 
is  due  solely  to  earth  rate. 

Figure  B-l  presents  the  F matrix  which  describes  the 
dynamics  of  all  the  states  appearing  in  the  error  analysis 
computer  program.  Table  B-2  defines  the  parameters  which 
appear  in  the  F matrix.  To  complete  the  specification  of 
those  matrices  necessary  to  propagate  the  covariance  matrix 
in  the  error  analysis  program,  it  is  necessary  to  define 
the  continuous  noise  matrix  Q that  appears  in  Eq.  (3.3-5). 
The  only  elements  of  the  Q matrix  which  are  non-zero  are 
the  diagonal  elements  corresponding  to  the  random  drift  rate 
and  accelerometer  errors.  For  these  states,  which  are  mod- 
eled as  first  order  Markov  processes,  the  value  entered  in 
the  Q matrix  is  2o  B where  1/8  is  the  correlation  time  of 
the  process  and  o is  the  rms  value  of  the  error. 

Because  many  of  the  states  in  the  MICRON  system  er- 
ror vector  model  the  equivalent  effects  in  the  L-frame  of 
errors  which  actually  occur  in  the  gyro  case  fixed  M-frame, 
the  initialization  of  the  covariance  matrix,  P,  defined  as 


fc*T] 


(B-l) 


is  not  as  simple  as  for  classical  inertial  system  error  analy- 
ses. Typically,  the  system  error  sources  are  assumed  to  be 
statistically  independent  resulting  in  an  inertial  covari- 
ance matrix  that  is  diagonal.  In  the  MICRON  system  analysis, 
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Figure  B-l 


0 

0 

Q 

b 

o 

0 

-% 

0 

"k 

“d 

0 

b 

b 

ft 

D 

“b 

o 

0 

*IflD 

0 

b 

an* 

b 

o 

0 

-s. 

0 

B* 

0 

b 

0 

0 

D 

igure  B-l 


F-Matrix  for  Error  Analysis  Program 


jrZLT’ 


TABLE  B-2  (Continued) 
PARAMETERS  USED  IN  FIGURE  B-l 


n - earth' 8 angular  rate 
- Q cob  (Lat) 

0D 0 sin  (Lat) 

g — gravity 


p — axial  mass  unbalance  drift  rate 
® compensation  coefficient 

R - radius  of  the  earth 


1/Bya*  1/®yE*  1^®VD 
I/&eN»  I/®eE* 


— vertical  channel  damping  coefficients 
(see  Appendix  A) 

— correlation  time  of  random  accelerometer 
errors 

— correlation  time  of  random  drift 
rate  errors 
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unit  vectors  along  case-fixed  X.Y.Z 
axes  coordinatized  in  the  L-frame 


for  those  error  sources  which  occur  in  the  M-frame,  the  as- 
sumption of  statistical  independence  was  made  relative  to 
this  frame.  Consequently,  there  was  considerable  cross- 
correlation among  the  corresponding  model  states  which  were 
defined  relative  to  the  L-frame.  Also,  because  many  of  the 
model  states  were  products  of  an  error  source  and  a function 
of  the  MESG  RSA  direction  vectors,  it  was  necessary  to  con- 
struct the  initial  covariance  matrix,  Pq,  from  the  program 
input  parameters  which  described  the  error  statistics  and 
the  RSA  vectors. 

To  compute  the  L-frame  covariance  from  the  M-frame 
statistics  requires  a simple  transformation  using  CM  , the 
coordinate  transformation  from  the  M-frame  to  the  L-frame. 
For  example,  if  Pe^  is  the  covariance  matrix  defining  the 
three  case-fixed  bias  terms  of  MESG  No.  1 in  the  gyro  (M) 
frame, 
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can  be  computed  as 
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PeBl  CM 


(B-4 ) 


The  process  of  constructing  an  initial  covariance 
matrix  corresponding  to  the  state  vector  defined  in  Table  B-l 
from  program  input  parameters  is  straightforward,  but  since 
it  is  rather  tedious  to  enumerate,  a detailed  presentation 
will  not  be  included  in  the  appendix.  With  the  definition 
of  the  system  state  vector  and  the  system  matrices  F,  Q and 
P,  all  the  equations  necessary  for  the  covariance  error 
analysis  program  have  been  presented. 
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appendix  c 


REDUNDANT  AXIS  CONTROL  IN  THE 
MICRON  MECHANIZATION 


The  purpose  of  this  appendix  is  to  demonstrate  that 
the  redundant  axis  control  that  appears  in  the  MICRON  mech- 
anization equations  (Ref.  14)  maintains 

-N1  ‘ — N2  = Icl  * Ic2  (C-l) 

where  §N1  and  SN2  are  the  estimates  of  the  RSA  direction 
vector  in  the  navigation  frame  coordinates  and  and  y 
are  the  compensated  readouts  of  the  RSA  vectors  (coordin^- 
tized  in  the  M-frame).  As  was  discussed  in  Chapter  3,  it 
is  important  that  Eq.  (C-l)  be  valid  because  the  inertial 
(S)  frame  in  which  the  accelerometer  outputs  are  integrated 
is  not  an  orthogonal  coordinate  system.  For  this  presenta- 
tion, the  updates  to  the  spin  vector  estimates,  SN1  and  S 
to  incorporate  the  effects  of  vehicle  motion,  earth  rateT^' 
and  estimated  drift  rate  are  ignored  since  they  are  applied 
correctly  to  each  direction  cosine  vector  and  are  thus  not 
an  important  part  of  the  redundant  axis  control  analysis. 

Considering  S»  as  the  primary  MESG , Fig.  C-l  is  a 
block  diagram  representation  of  the  redundant  axis  control. 
Ignoring  the  sampled  data  aspects,  since  the  1/16  sec  sampl- 
ing rate  should  not  introduce  any  significant  error,  it  can 
be  seen  that 
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Figure  C-l  Block  Diagram  of  Redundant  Axis 
Control  (MESG  No.  1 Primary) 


5ni  * — N2 


— N1  * (^N2  + A-N2) 

— N1  * ( — N2  + T*k  (*cl  * *c2  " — N1  * 2 ) ^Nl) 

— N1  * — N2  “ I7ic(-Nl  * ^2) + I+k(Ic1  * ^2) 

(C-2) 


Equation  (C-2)  can  be  seen  to  force  Sjq  * to  track 

. y 2 which  is  the  desired  effect  of  redundant  axis  con- 
. The  MICRON  mechanization  has  selected  the  gain 


1 - (Sni  ' Sn2)2 


(C-3) 
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The  choice  oi  gain  in  Eq.  (C-3)  increases  the  redundant  axis 

centre!  in  a manner  proportional  to  the  nonorthogonaity 

S and  S As  the  nonorthogonality  increases,  it  is  more 

critical  that  the  estimated  angle  between  the  spin  vectors 

re  • S ) track  the  measured  angle  (Ycl  * Ic2^  t0  aV° 

'“-Nl  -N2  f the  n0northogonal  in- 

introducing  errors  due  to  the  use  oi 

ertial  (S)  coordinate  frame  in  the  MICRON  mechanization. 
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